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Introduction 



This thesis is mainly about how to set up and carry out the idea of back-reaction in 
a physically meaningful way. 

Most, if not all, observations are consistent with the cosmic concordance model 
according to which, today, one-fourth of the mass-energy of the universe is clustered 
and dominated by cold dark matter. The remaining three-quarters is uniform and 
dominated by a fluid with a negative pressure (dark energy, or A). 

While the standard ACDM model seems capable of accounting for the observa- 
tions, it does have the feature that approximately 95% of the mass-energy of the 
present universe is unknown. We are either presented with the opportunity of dis- 
covering the nature of dark matter and dark energy, or nature might be different 
than described by the ACDM model. Regardless, until such time as dark mat- 
ter and dark energy are completely understood, it is useful to look for alternative 
cosmological models that fit the data. 

One non-standard possibility is that there are large effects on the observed ex- 
pansion rate (and hence on other observables) due to the back-reaction of inho- 
mogeneities in the universe. The basic idea is that all evidence for dark energy 
comes from the observational determinations of the expansion history of the uni- 
verse. Anything that affects the observed expansion history of the universe alters 
the determination of the parameters of dark energy; in the extreme it may remove 
the need for dark energy. 

The "safe" consequence of the success of the concordance model is that the 
isotropic and homogeneous ACDM model is a good phenomenological fit to the real 
inhomogeneous universe. And this is, in some sense, a verification of the cosmological 
principle: the inhomogeneous universe can be described by means of an isotropic and 
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INTRODUCTION 



homogeneous solution. However, this does not imply that a primary source of dark 
energy exists, but only that it exists as far as the phenomenological fit is concerned. 
For example, it is not straightforward that the universe is accelerating. If dark energy 
does not exist at a fundamental level, its presence in the concordance model would 
tell us that the pure-matter inhomogeneous model has been renormalized, from the 
phenomenological point of view (luminosity-distance and redshift of photons), into 
a homogeneous ACDM model. 

There are two ways to approach the problem of the back-reaction. 
The first point of view is strictly observational and focuses directly on the past light 
cone, on the effects of large-scale nonlinear inhomogeneities on observables such as 
the luminosity-distance-redshift relation. 

The second point of view tries, instead, to interpret the inhomogeneous universe by 
means of an effective model. 

These two different points of view actually share the same idea of smoothing out 
the inhomogeneities. The differences are in the way this process is carried out. 
In the approach which uses an effective model the averaging is explicit, while in 
the observational approach the averaging is implicit: inhomogeneities are indeed 
"integrated out" in the luminosity-distance-redshift relation. 

This duality will shape the structure of the present thesis: we will present the 
theoretical backgrounds and work out illustrative models for both the approaches. 

After an outline of the evidences for dark energy. Chapter [T] will introduce the 
candidates we are interested in: the cosmo logical constant, the quintessence and the 
back-reaction. This will give us the opportunity to give our point of view about 
the dark energy problem and to justify the idea of the back-reaction as a possible 
alternative explanation for dark energy. 

Chapter [2] will be about the dual point of view we just talked about. In line with 
Ellis's work, we will first introduce the smoothing process in general relativistic 
cosmology and then we will try to work it out theoretically by setting the problem 
in a physically meaningful way. 

In Chapters [3] and H] we will attempt to work out the theoretical concepts in- 
troduced by means of concrete models, even if it is difficult to analyze a realistic 
inhomogeneous model of the universe in an exact way. In order to understand the 
physics of the back-reaction, therefore, we have looked at exact general relativistic 
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solutions where some symmetries will enable us to carry out the calculations. We 
have chosen the spherically symmetric dust Lemaitre-Tolman-Bondi (LTB) solution 
as our starting point. 

There is increasing attention to the LTB model because it was shown that it is 
possible to fit the observed luminosity-distance-redshift relation by adjusting the 
LTB free functions. To achieve this result, however, it is necessary to place the 
observer at the center of a rather big underdensity. Even though we will build 
nothing more than toy models, we will try to avoid this fine-tuning. The ultimate 
goal is indeed to find a realistic dust model that can explain observations (like the 
luminosity-distance-redshift relation) without a need for dark energy. The ultimate 
aim being to have an exactly solvable realistic inhomogeneous model. The models 
we will build are but a first small step in this pragmatic and necessary direction. 

In Chapter [3] we will study the average scheme developed by Buchert building 
a "homogeneous" universe model where the parameters are chosen in order not to 
single out the center. We will study both the curved and the fiat case. 

Then, in Chapter HI we will build a Swiss-cheese model with the observer in the 
cheese looking through a series of holes. The cheese consists of the usual Friedmann- 
Robertson- Walker (FRW) solution and the holes are constructed out of a LTB so- 
lution. We will study this model under both points of view. First, we will focus 
on the effects of large-scale non-linear inhomogeneities on observables such as the 
luminosity-distance-redshift relation. Then we will try to apply the fitting scheme 
developed by Ellis and Stoeger: we propose a fitting procedure that is intermediate 
between the fitting approach and the averaging one: a fit with respect to light-cone 
averages. 

In Appendix |X] we will present a general model for the cosmological evolution of 
the fine structure constant driven by a typical Quintessence scenario. We consider 
a coupling, between the Quintessence scalar and the electromagnetic kinetic term, 
given by a general function. We study the dependence of the cosmological variation 
of the fine structure constant upon the functional form of the general function chosen 
and discuss the constraints imposed by the data. We find that different cosmological 
histories for the fine structure constant are possible within the available constraints. 
Appendices |B] and O are, instead, about LTB models. 

For a discussion of the results see the Conclusions. 
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Chapter 1 
Dark Energy 



In this chapter we will introduce the Dark-Energy problem, at first its experimental 
evidences and then the possible explanations in which we are interested. It will help 
us to introduce the Back-Reaction problem in the next chapter. 

1.1 Evidences 

The fact that matter energy density does not dominate the universe has been one 
of the most stunning discoveries in the last few years in cosmology. 

Most, if not all, observations are consistent with the cosmic concordance model 
according to which, today, one-fourth of the mass-energy of the universe is clustered 
and dominated by cold dark matter. The remaining three-quarters is uniform and 
dominated by a fiuid with a negative pressure which we call dark energy. 

The three main experimental evidences for dark energy are summarized in Fig. ll.li 
The blue area is about the results of the 2dF galaxy survey which gives — 0.3. 
This measure is independent of fir, e and so it is vertical^] in the plane flM-f^DE- 
The orange area represents the constraints from the WMAP observations of the 
CMB anisotropics. They give 1 ^ flror = flhi + floE and so the orange area is 
aligned to the line for (0, 1) and (1, 0). 



^The blue area is not actually vertical, but slightly tilted. This because the 2dF galaxy survey 
is mainly at z = 0.14. In order to evolve the data till z = a cosmological model (in particular a 
ACDM model) has been used and therefore a dependence on dark energy introduced. 
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2 CHAPTER 1. DARK ENERGY 

Finally, the green area shows the Supernovae la measurements about the de- 
celeration parameter which, in the ACDM model, is given by q = VLm/'^ — ^de- 
Therefore the green area is almost perpendicular to the orange one. 

The fact that these constrains cross perpendicularly and consistently gave the 
name cosmic concordance model to a ACDM with: 

nM ^ 0.25 
^ 0.75 

wde ^ -1 (1-1) 




Figure 1.1: Experimental constraints on the plane Qm-^de from the 2dF galaxy survey 
(blue), CMB (orange) and Supernovae la (green). 
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1.2 Candidates 



In this section we will take into account three possible explanations for dark energy: 
the cosmological constant, a cosmological scalar field and inhomogeneous universe 
models. In this thesis we are really interested only in the last possibility. We will 
describe the first two because, beside being interesting in themselves, they will help 
us understand the problem and will turn out to be useful to describe inhomoge- 
neous models effectively. They are indeed based on the cosmological principle of 
homogeneity and isotropy. 

For a more comprehensive examination of possible approaches and different 
points of view to dark energy see [19]. 

1.2.1 The Cosmological Constant 

Trying to explain dark energy through the cosmological constant A, that is by us- 
ing the concordance model, has the advantage of simplicity and agreement with 
experimental data. 

The cosmological constant appears as a free parameter in the Einstein's equa- 
tions: 



where we are using the signature (— , +, +, +) and geometric units, c = 1 = G. 

The freedom about the value of A is cause of problems. In spite of being treated 
as the single problem of the value of A, from a conceptual point of view they are 
actually three distinct problem^. This distinction will help us understand which 
ones are more pertinent to the cosmological problem and which are actually the 
advantages in using alternative explanations for dark energy. 

For a wider overview of cosmology in the presence of the cosmological constant 
we refer to [65l [HI EDI ESI El] . 

Initial conditions problem 

The first problem is about the initial conditions that we have to give to the con- 
cordance model: it would be expectable to have comparable values of the initial 

^For a similar approach see [57] , 
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Figure 1.2: Evolution of the energy density of radiation (short-dashed hne), matter (long- 
dashed line) and cosmological constant (solid line) with respect to the redshift. All energy 
densities are expressed in units of the present critical energy density pco ■ 

energy densities. However, as you can see from Fig. 11.21 matter and radiation ini- 
tial energy-density values are much bigger than the cosmological constant one. The 
concordance model implies therefore a sizeable fine tuning of the initial conditions. 

Coincidence problem 

The second issue is about the values of the present-day energy density of matter 
and cosmological constant. As you can see from Fig. 11.21 they evolve differently 
and therefore generally, during the evolution of the universe, we will expect them 
to have different values. However their values happen to be of the same order at 
present time. This fact needs a sizeable fine tuning to happen. 

Origin problem 

Finally, problems show up if we ask where the cosmological constant comes from. 

Classically, it is a free parameter given as an initial condition, p^- 

Quantum mechanically, we have the contribution from the ground state of the 
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energy- momentum tensor which does not need to be zero: 

{0\T,,\0)^T;:f = PGsg,u (1.3) 

Pes has the same equation of state of p\ and therefore we can combine them into 
the vacuum energy density pv- 

PV = PA + PGS (1-4) 

The problem is that any quantum-mechanic estimation of pv largely exceeds the 
experimental value of pde ~ (10^^ eV)^. Typical estimations are around pv ~ 
(10^'' eV)^, that is, 30 orders of magnitude bigger. 



1.2.2 Quintessence models 

Quintessence or cosmological scalar-field models will be, as effective models, a useful 
tool in our analysis of a inhomogeneous universe. 

The idea is to let the cosmological constant be a dynamical quantity. The La- 
grangian of a quintessence field is, therefore: 



£^ = -9^09^0-^(0) (1.5) 



from which, assuming spatial homogeneity, we obtain: 



^02 + y(0) e p^ = l02-T/(0) (1.6) 



The equation of state is: 



Pi, 4>y2 + v{4,) 

Quintessence models are extensively studied (see for example [26]) mainly because of 
their flexibility. In particular there exist classes of potentials which exhibit attractor 
solutions independent of initial conditions [6pj. 

For illustrative purposes, in Fig. I1.3l we show the evolution of the energy density 
and of the equation of state for a scalar potential V = M^(f)~^ and initial conditions 
PjT/Pc'o = 10'^° at z = 10^°. M gives the energy scale of the potential and has to be 
fine-tuned in order to have the right present-day quintessence energy density. 
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As you can see from Fig. I1.3[ at early times the kinetic energy is dominant 
over the potential energy: during this period, called kination and characterized by 
W0 ~ 1, the scalar field rapidly rolls down along the potential. Then the friction 
performed by the Hubble terno slows down the run of the scalar till its freeze: during 
this period, called slow-roll and characterized by w,^ ~ —1, the potential energy 
dominates and the scalar mimics a cosmological constant. Finally the scalar reaches 
the attractor and follows a dynamics characterized by its features. In the illustrative 
case shown, the attractor is characterized by ~ —0.67. See Appendix El for more 
details. 




Figure 1.3: Evolution of the energy densities (left) and scalar equation of state (right) for 
a quintessence model with potential V = M^cf)^^ and initial conditions p^^ / pco = 10'^'^ at 
z = W^^. The dot-dashed line represents the energy density of radiation, the dotted line 
the energy density of matter, the green dashed line the energy density of quintessence and 
the red solid line the attractor. All of the energy densities are expressed in units of the 
present critical energy density pco- From J45|. 



Pros 

A first advantage of the example shown is that it solves the initial conditions prob- 
lem: we can change of many order of magnitude the chosen value of p^"/ pco = 10^° 

■^The continuity equation for a quintessence field can be written as — —'SH(p^. 



1.2. CANDIDATES 



7 



at z = 10^° and still the late time dynamics will be the same: only the moment at 
which the scalar reaches the attractor will change. This is a real improvement in 
comparison with the cosmological constant. 

In the second place, a dynamical scalar field may have chances to solve the co- 
incidence problem. Indeed, if it couples to the matter energy density, it could be 
possible to explain why, today, they are of the same order of magnitude. However, 
quintessence models generally feature a sizeable fine-tuning about the mass param- 
eter of the potential, similarly to what happens with the cosmological constant. 

Cons 

Beside the problem of the origin of quintessence, its first drawback is about its mass 
which has to be ultra-light: 

m^^ ~ ifo ~ 10"^^ eV (1.8) 

Quantum-mechanics direct or indirect corrections will easily bring m^^ toward much 
bigger values. 

Other problems come from the possible coupling of the quintessence field with the 
other terms in the matter-radiation Lagrangian and from the fact that such an ultra- 
light field has long-range interactions: A ~ m'^^ ~ H^^. Variations of fundamental 
constants and violations of the equivalence principle are therefore expected. For 
more details see Appendix [XI 

Effective theory? 

In the case quintessence models are used as effective models, none of the drawback 
examined above will be applicable. They indeed are about the quantum mechanics 
implications a primary quintessence field brings to light. 

As for this thesis, the most appealing feature of quintessence is its dynamical 
nature which can be linked to the evolution of the universe. In particular, a cou- 
pling of an effective quintessence field to inhomogeneities could solve the coincidence 
problem. We will discuss this point in the next chapter. 
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1.2.3 Inhomogeneous models - a point of view 

This overview of advantages and drawbacks of ACDM and quintessence models gives 
us some hints to better set the dark energy problem. We will present now the point 
of view that will characterize this thesis. It is clarifying to write: 

Pv = (pa +) Pgs ^ Pde ~ Pm (1-9) 

On the left-hand side, the two terms p\ and pos are really different and could be 
inappropriate to mix them. The latter comes from quantum mechanics, the former 
comes from general relativity. p\ is likely connected to cosmology and therefore to 
Pde while for pos it could not be the case. Indeed typical estimations for pas are 
far away from the cosmological value. 

Therefore it could be possible that we should be only concerned about pa the role 
of which has to be understood in order to have a satisfactory general-relativist 
cosmological model. 

On the right-hand side, the similarity in density values of the dark energy and 
matter components suggests a direct or indirect connection: the dark energy problem 
might be a cosmological one. 

If we accept this reasoning and remember that the homogeneous ACDM model 
is a good fit to the real inhomogeneous universe, we could think that the connection 
is by means of the smoothing of the inhomogeneities: the back-reaction of inhomo- 
geneities makes an inhomogeneous matter FEW model appear as a ACDM model. 
One immediate benefit is to turn the coincidence problem into a hint in favor of the 
back-reaction: only recently have the cosmic structures evolved enough to have a 
sufficiently strong back-reaction of inhomogeneities which could affect observables. 

The issue now is what we mean by back-reaction of inhomogeneities. There are 
broadly speaking two distinct approaches as we will now discuss. One is focused on 
theoretical aspects while the other on observations. 

Both the approaches will, however, share the idea of smoothing out inhomogeneities. 
The duality in the interpretation of this concept will characterize this thesis. 
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Observational side - Pde 

The "safe" consequence of the success of the concordance model is that the isotropic 
and homogeneous ACDM model is a good observationalfit to the real inhomogeneous 
universe. And this is, in some sense, a verification of the cosmological principle: the 
inhomogeneous universe can be described by means of an isotropic and homogeneous 
solution. 

However this does not imply that a primary source of dark energy exists, but only 
that it exists effectively as far as the observational fit is concerned. For example it is 
not straightforward that the universe, as Supernovae la measurements seem to tell 
us, is globally accelerating. If dark energy does not exist at a fundamental level, its 
evidence coming from the concordance model would tell us that the purely-matter 
inhomogeneous model has been renormalized, from the observational point of view 
(luminosity and redshift of photons), into a homogeneous ACDM model. Moreover, 
the very homogeneous nature of dark energy seems a clue about its effective nature. 

There could indeed be the possibility that there are large effects on the observed 
expansion rate due to the back-reaction of inhomogeneities in the universe (see, e.g., 
Ref. [36] and refs. therein). The basic idea is that all evidence for dark energy comes 
from observational determination of the expansion history of the universe. Anything 
that affects the observed expansion history of the universe alters the determination 
of the parameters of dark energy; in the extreme it may remove the need for dark 
energy. 

Summarizing, this approach is tied to our past light cone: it will focus on the ef- 
fects of large-scale non-linear inhomogeneities on observables such as the luminosity- 
distance-redshift relation. Even though there is no explicit averaging in here, the 
dj^iz) is a about the luminosity and redshift of photons that travelled through in- 
homogeneities: their effects are therefore "integrated out", averaged in an implicit 
way. 
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Theoretical side - pA 

We said that pA is likely connected to cosmology and that, in any case, we have to 
understand its meaning in order to have a satisfactory general-relativistic cosmolog- 
ical model. 

It could be that pA is not directly related to what we measure as pde- The latter 
comes from the fitting of the real inhomogeneous universe by means of a ACDM, 
while the former is about a free parameter of Einstein's equations. 

As we will discuss in the next chapter, pA could be the result of a smoothing 
process [311 ES]- This is an important issue in general relativity because, generally, 
every measurement involves some form of smoothing or averaging. 

Following this point of view, pA could be a product of the back-reaction of in- 
homogeneities and have its freedom fixed by a measuring process. Inhomogeneities 
introduce a scale in the otherwise scale-free general relativity. We think this is a 
crucial step in understanding how General Relativity effectively works in a lumpy 
universe. 

As we will see, the issue is to understand how to carry out this smoothing 
meaningfully. In particular if it is possible to connect pA to Pde- Within the previous 
approach we were only concern about the luminosity-distance-relation, while here 
we want, in addition to that, a phenomeno logical model that fits observations, in 
other words, we want a description by means of a mean field. 



Chapter 2 

The Back-Reaction Problem 



In this chapter, after setting the problem in the first section, we are going to explore 
the dual point of view sketched in Sect. I1.2.3I 

The second section will suggest two possible approaches to the smoothing pro- 
cess. The approach based on Buchert's equations will be carried out in Chapter El 
while the one based on the light-cone fitting will be studied in Chapter |H 

The third section will be, instead, about the strictly observational point of view 
which will be worked out in Chapter HI 



2.1 Foundations: the smoothing process 

We will here set the smoothing process. We think that this is a crucial step in 
understanding how General Relativity effectively works in a lumpy universe. 

This section will lay the foundations of the problem without actually telling how 
to pragmatically search for a back-reaction effect. The following two sections will 
be about finding out how to achieve the smoothing process. To set the problem we 
will follow [3III32]. 
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CHAPTER 2. THE BACK-REACTION PROBLEM 



Let's look at fig. 12.11 as a starting point: it compares models of the same re- 
gion of the universe on three scales showing different amounts of detail. Scale 1 
represents all details down to stars. Scale 3 represents all details down to galaxies. 
Scale 5 represents large scale features only. The usual models of cosmology, like the 
concordance model, are shown at Scale 5. It should be pointed out that these dif- 
ferent matter tensors and metric tensors are intended to describe the same physical 
system and the same space-time, but at different scales of description. Smoothing 
out inhomogeneities renormalizes the description of the universe. 



Density 



A |SCALE1| (M,,g,,,,T,,, 



distance 



S31 

J. smoothing 
^ ' map 



Density 
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^ ' map 



distance 



Figure 2.1: Comparison of models of the same region of universe on three scales showing 
different amounts of detail. Scale 1 represents all details down to stars. Scale 3 represents 
all details down to galaxies. Scale 5 represents only large scale features. From j3Jf . 



2.1. FOUNDATIONS: THE SMOOTHING PROCESS 



13 



General relativity tests confirm that Einstein's equations hold on Scale 1 which 
is the starting point in the fiow chart below. 



9l ab 



ab ~ Mqb 
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The fiow chart summarizes tensors correspondence from Scale 1 to Scale 3. The 
two models are related by: 

• Maps S*, which determines which points in the different underlying manifolds 
are related to each other by the smoothing procedure. 

• Maps S, which determines the metric tensor of the smoothed-out representa- 
tion from the more detailed one. 

• Maps 5", which determines the matter tensor. 

• Maps S", which determines the Einstein tensor. 

Because of the non-linear nature of the fields equations, in general, the operations 
of smoothing will not commute with going to the field equations, that is: 

S" ^ S' or {G,,{g^^)) + (2.1) 

We can better set the problem defining a tensor P3 representing the difference be- 
tween the Einstein tensor 6*3 defined from the smoothed-out metric g^a, and the 
smoothed-out matter tensor T-^. This correction will take care of the change of scale 
of description. 

GzaP = RsaP — 2 93 aP -^3 = SttTsq,^ + SttPsq/J (2.2) 

The tensor P3 represents the effects of small-scale inhomogeneities in the universe 
on the dynamic behavior at the smoothed-out scale. 
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Comments 

The implicit assumption in the usual approach is that at Scale 5 the term P5 vanishes. 
Clearly this may not be true. Moreover, P does not need to obey the usual energy 
conditions. It is indeed an effective term. There would be no conceptual problems 
in having wp < —1, for example. 

It is difficult to refrain from considering the possibility of a connection between 
the cosmological constant and this correction term: 

Pal3 = PAgal3 that is Wp = -1 (2.3) 

Pa could be indeed a product of the back-reaction of inhomogeneities and have its 
freedom fixed by the smoothing process which, generally, is introduced by every 
measurements. Inhomogeneities introduce a scale in the otherwise scale-free general 
relativity. 

Even if equation (12.31) does not exactly hold, it still suggests the idea that the 
correction term P may be the source for an effective quintessence model. In the 
next section we will show a concrete example about this. 
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2.2 Smoothing out inhomogeneities - theoretical 
side 

In the last section we set the foundation of the problem of the back-reaction. Now 
we will consider the theoretical approach in smoothing out inhomogeneities, that 
is, we will look for a phenomenological model that fits observations, that is, for a 
description by means of a mean field. 

The two approaches we are going to describe will give two conceptually different 
PA-like terms. The problem is really to set the back-reaction problem in a physically 
meaningful way. 

2.2.1 Averaging the traces of Einstein's equations 

We will describe now a way to average Einstein's equations based on averaging their 
traces 

Equations 

The matter content be made of irrotational dust with four-velocity = (1, 0, 0, 0), 
that is, the local observer is comoving with the energy flow of the fluid. Hence 
T^u = P u^Uu and the Einstein equations are: 

Rfiu - ^9fiuR = 8npUf,u^ (2.4) 

where p is the matter energy density and R the Ricci scalar. We can work in the 
synchronous and comoving gauge with line element: 

ds"^ = -dt^ + Qij dx'dx^ (2.5) 

where t is cosmic time. Greek indices run through 0...3, while Latin indices run 
through 1...3. 

We will introduce the expansion tensoi0 9*^ : 

Qij = h"^ h^- Ua-p = Ui-j where /i"^ = + (2.6) 
^The extrinsic curvature is the opposite 
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from which we can define the volume expansion scalar 9 and the traceless tensor 
cr*^, called shear: 



- + (2.7) 



9 
3 

Now let us define the averaging operation which for a scalar field is a covariant 
operation, given a foliation of spacetime: 

{:F)n = ^ [:FdV= Yf,t ^2.8) 

where the average depends on content, shape and position of the comoving spatial 
domain of averaging, which is considered given. 

We will introduce a dimensionless effective scale factor via the volume: 

Thus, the averaged expansion rate may be written in terms of the scale factoi@: 

{9)d = ^ = 3^^3Ho (2.10) 
Vd cld 

Dots denote derivatives with respect to cosmic time. 

The smoothing procedure and the dynamical evolution does not commute as it 
is clear from this relation: 

(^)l, - {J')d = {:F9)d - {9)d{:F)d (2.11) 

Within this averaging scheme, this is the source of the backreaction effect. Indeed, 
for JF = 6* we will have: 

WL = (^)D + (^%-(^)i> (2-12) 

and therefore even if {9)jj < 0, it may happen that (9)'^ > 0, that is the inhomo- 
geneities allow the coarse-grained deceleration parameter 

_ 3(g)p + {9)1 _ io 
to be negative in spite of its positive local valu^. 



^It follows from J — 6 J where J — ^/g 

^It follows from the local Raychaudhuri's equatfon (|2?T5l) that q = -{39 + 6'^)/0'^ = 6{<j^ + 
2iip)/e'^ > and that {e)jj < 



2.2. SMOOTHING OUT INHOMOGENEITIES - THEORETICAL SIDE 



17 



Averaging 

Now we can average the following local scalar equations: 

II + -7^ — Stt p = energy constraint (2-14) 
I + I^ — 2 II + 47rp = Ray chaudhuri's equation (2-15) 
p + Ip = continuity equation (2.16) 

where two of the scalar invariants of the expansion tensor, namely its trace and the 
dispersion of its diagonal components, are: 

I = e', = e 

II ^ l{e'-e',e\) = le^-a' (2.17) 

where = ^a'^^a-'^. The IZ in fl2.14l) . and in the following equations, is the trace 
of the spatial Ricci tensoi0 TZij. Actually, the density is a scalar invariant as well, 
T = g^'^T^u = —p- The result of the averaging is: 

(II)B + ^(7^)^-87r(p)^ = (2.18) 
(1)1, + (I)i- 2(11)^ + 47r(p)^ = (2.19) 
{p)h + {i)D{p)D = (2.20) 

that is, provided we express the equations in terms of the invariants fl2.17p . the 
averaged quantities obey the same equations as the local ones in spite of the non- 
commutativity of the averaging procedure and the dynamical evolution, which is 
expressed by the commutation rule (12.111) . 

The reason for this nontrivial property is the special type of nonlinearities featured 
by the gravitational system, in particular the nonlinearity in 6 contained in Ray- 
chaudhuri's equation. 

Summarizing, the traces of Einstein's equations, expressed by means of the scalar 
invariants I, II and p, are not structurally affected by the averaging procedure: the 
corrections are in the renormalized invariants and not in their relations. This result 
has some consonance with the idea of pa fixed by the smoothing process: not even 
in that case Einstein's equations are structurally changed. 

TZij is the Ricci tensor of the metric gij of the hypersurface and not Rij — /i" h^j Rap- 
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The extra-terms produced by the averaging will appear if we rewrite the first 
two equations in a more familiar form: 

Hi = y(p)D-^(7^)^-^QD (2.21) 

cld 47r , , 1 ^ , , 

— = -^{p)d + -Qd (2.22) 

where the (kinematical) back-reaction source term is: 

Qn ^ 2{U)n - 1(1)1 = ^ - {6)1) - 2{a')n (2.23) 

The equations fl2.2m2.2^ can once more be rewritten in a dimensionless form: 

1 = fif, + r]^ + fig (2.24) 

qn = ^n^t + ^^S (2-25) 

where we have defined the density parameters Q'^j = 87t{p) j-^/SHl, Qq = —Qd/GHI 
and = — (7^)^/6//|, in the usual way. 

The consistency of equations (12. 21112.221) requires that Qn and (T^)^) satisfy the 
integrability condition: 

{a%Qn)' +a%{al{n)n)' = (2.26) 
Effective back-reaction source 

Finally we can rewrite again (12.21112.221) in term of the effective backreaction source: 

Hi = Y((p)z, + PBi?) ^ ^1^, + ^^^ = -^ (2.27) 

^ = + (2.28) 

an o '—^ 

% 

where the effective backreaction source energy density and pressure are: 

lOTT lOTT lOTT O lOTT 

The equation of state will be: 

_ Vbr _ -Qd + (7^)^/3 

- - -Qn - in), ^'-''^ 
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We have therefore the non trivial result that the standard Friedmann equations are 
modified only by an extra source beside their dependence on the scale of averaging. 
Unfortunately the system of the equations (12.27112.2^ is not closed, the smoothing 
procedure has washed out too many details. To close the system it will be sufficient 
to have the equation of state: we will come back to this point later in this section. 

These last equations look like the corresponding ones for quintessence, f ll.5til.6l) . 
with kinetic energy ~ —Qd and potential ~ —(7^)^. For example, with zero 
potential (i.e. (7?.)^ = 0) we have from (12.261) ~ 1/ a]f which is the usual kination 
with w = 1. Indeed a topic of research could be to analyze the problem with the 
quintessence formalism, trying to find potentials coupling (p to the inhomogeneities 
or some other connected observables. Because (p is an effective field, it does not 
suffer of the problems examined in Sect. I1.2.2[ It will be possible also to have 
w < —1. See for a work in this direction. 

Acceleration in a dust universe 

The condition to have acceleration follows from the acceleration equation (12.221) and 
(ESS}: 

-fig>^ or Qd>^7i{p)d (2.31) 

However, having significant backreaction is not enough: it is necessary that 
and (TZ)^ are strongly coupled, otherwise Qd will fall in kination and will become 
quickly negligible. Therefore (11)^ will evolve in a different way with respect to the 
usual curvature k: only if Qd and {71)^ are decoupled it will be {71)^ ~ V'^d- So 
the CMB tells us only that the part of (7?.)^ which evolves as l/a|) is negligible. 

Assuming = 0.3 we have plotted in fig. I2.2l the value of the main parameters 
as a function of Qq to see if they are in contrast with the present day data. From 
the plot we see that the following values: 

- 0.45 < fig < -0.3 and 1 < fig < 1.2 (2.32) 

are in agreement with the data. We stressed again that the value we have found for 
fig is not ruled out by the CMB because the latter constrains only the part of fig 
which evolves as l/a|). The volume D refers to the scale at which the parameters 
are measured. 
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1 . 4 
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Figure 2.2: We have plotted for = 0.3 the value of the principal parameters as a 
function of Qq. The dashed line is i^^, the dotted line is qo SLnd the dot-dashed line is 
wbr- The measured values J59| are marked in solid red, while in the gray area there is no 
acceleration. 



A conjecture 

We can try to apply the conjecture that = pA da/B to this formahsm, that is: 

WBR = -1 (2.33) 
We are now able to close the system (12. 27112. 28|) : 

~ ~Q^^ ~ const(a£)) or Qd = —i:{T^)d ~ const(aD) (2.34) 

o o 

where the last equalities derive from the Constance of pbr which is implied by (12.331) . 
The solution (I2.34p satisfies the integrability condition (12.261) . 

It is possible to put a remedy to our ignorance about the dependence of the 
observable from qd expressing the scale as a function of the matter parameter: 
dD = ciD{^M)y so that to give is as to give the domain of averaging. We can 
then plot in fig. 12.31 this exact solution: 

^S = l^M-l ^n=l-l^M <lD = lnZ-l (2.35) 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 2.3: Assuming wbr = —1, we have plotted the principal parameters as a function 
of J^^/. The dashed line is ft^, the solid line is 0,q and the dotted line is qd- 



Summarizing: 



= 0.3 l^go = 0.7 Wbr = -1 Qd = -0.55 



= 0.3 


''''BR 


= 0.7 


where 


0^ 


= n^ + n 


and 




= -0.35 



(2.36) 



= 1.05 



We are stressing again that the volume D refers to the scale at which the parameters 
are measured. 



Bare and dressed cosmological parameters 

There is one more remark about the theoretical framework just described: the pa- 
rameters of eq. fl2.24p . which is: 

^^D^^D^-^ (2.37) 

are not directly accessible to observations because the analogous FRW ones implic- 
itly assume a constant curvature space [T6j . 
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Before making this comparison we need to average the inhomogeneous geometry 
into a constant curvature one. The following sketch gives an idea of the process. It 
seems as if we have to correct only for volume effects when dressing the parameters. 




However, averaging the curvature, we will find: 

/V-\ 

(7^)5 =(^j {'^)d-Qd (2.38) 

that is, beside volume effects, there is an extra correction named curvature back- 
reaction Q^. Its analytical expression is similar to the kinematical back-reaction Qd 
but there are, instead of the scalar invariants of the extrinsic curvature, the ones of 
the intrinsic curvature. 

From eq. fl2.38p we can understand the physical content of geometrical averaging. 
It makes transparent that, in the smoothed model, the averaged scalar curvature is 
'dressed' both by the volume effect and by the curvature backreaction effect itself. 
The volume effect is expected precisely in the form occurring in fl2.38p . if we think 
of comparing two regions of distinct volumes, but with the same matter content, in 
a constant curvature space (remember that a constant curvature space is propor- 
tional to the inverse square of the radius of curvature, hence the volume-exponent 
2/3). Whereas the backreaction term encodes the deviation of the averaged scalar 
curvature from a constant curvature model, e.g., a FLRW space section. 

Therefore the relation we have to use to compare results to experimental data 

is: 

+ + Og = 1 (2.39) 
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2.2.2 The fitting problem 

The other theoretical approach in smoothing out inhomogeneities we will consider 
is based on the fitting problem [32] . 

There are broadly speaking, two distinct approaches, which have been applied 
to understand the large-scale structure of the universe. 



(a) 




Figure 2.4: (a) An exactly uniform and spherically symmetrical FRW universe U' mapped 
into the lumpy universe U so as to give the best possible ht. (b) An exactly spherical sphere 
fitted to the lumpy world to give the best ht possible. From \3^ . 

The standard approach is to make the assumption of spatial homogeneity and 
isotropy on a large enough scale. Then it follows that the universe is represented by 
a FRW model. 

The main problem of this approach is that it simplifies the way the real lumpy uni- 
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verse should be averaged. It does not deal with the correction term P introduced 
in the previous section. 

The concordance model, however, fits experimental data very well: the direct con- 
sequence of its success is, indeed, that the isotropic and homogeneous ACDM model 
is a good observational fit to the real inhomogeneous universe. And this is, in some 
sense, a verification of the cosmological principle of of spatial homogeneity and 
isotropy on a large enough scale: the inhomogeneous universe can be described by 
means of a isotropic and homogeneous solution. However this does not imply that 
a primary dark energy component really exists, but only that it exists effectively as 
far as the observational fit is concerned. 

If dark energy does not primary exist, its evidence coming from the concordance 
model would tell us that the purely-matter inhomogeneous model has been renor- 
malized, from the observational point of view (luminosity and redshift of photons), 
into a homogeneous ACDM model. 

So, following this reasoning, it is not true, as stated in the discussion about the 
previous section, that the standard approach assumes P to vanish. The standard 
approach simply takes the P correction term as a primary source in the energy- 
momentum tensor. 

The other approach is to make no a priori assumptions of global symmetry and 
build up our universe model only on the basis of astronomical observations. The 
main problem with such an approach is the practical difficulty in implementing it. 

An approach which is intermediate between the two outlined above is based on 
the fitting procedure. It asks the question about which is the FRW model which 
best fits our lumpy universe. Fig. 12.41 (a). This question will lead to a procedure 
that will make us better understand how to interpret the large-scale FRW solution. 
The approach resembles that used in geodesy, where a perfect sphere is fitted to the 
"pear-shaped" earth; deviations of the real earth from the idealised model can then 
be measured and characterised. Fig. 12.41 (b). 

The best-fit will be implemented along the past light cone. This because a mean- 
ingful fitting procedure should be related directly to astronomical observations. 
The FRW model we have in mind to use in the fitting is a quintessence-like model 
with an effective source with varying equation of state. We will develop these ideas 
in Section ITUl 
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Let's start discussing the approach about averaging the traces of Einstein's equa- 
tions. This scenario is theoretically consistent, but there are two main issues. 

Is observationally relevant? 

We have to make it clear if this approach gives an observationally meaningful back- 
rection. We have to understand how the Pa ^ Pbr that comes from this scenario is 
related to what we measure as Pde- 

A first important observation is that we should be mainly concerned about the 
luminosity-distance-redshift relation dL{z) and look at the effects inhomogeneities 
have on it. Therefore it seems approximate to average at constant cosmic time, we 
should instead average along the light cone. 

To overcome this problem we will use the approach sketched in Sect. 12.2.21 As we 
will present in Chapter HI we will use a fitting scheme based on light-cone averages. 

The scenario developed by Buchert, beside giving precious hints toward the un- 
derstanding of the back-reaction problem, can, however, be conceptually better un- 
derstood if we regard the spatial average as an ensemble average [36]. 




The inhomogeneities are seen as variables that take random values over different 
realizations of a volume D at a fixed location or over different locations of the vol- 
ume D in a single realization (fair sample hypothesis). In order to proper represent 
an observer in the probability ensemble or a randomly placed observer, we should 
average his measurements over the possible realizations. Averaging at constant cos- 
mic time over a volume of size comparable with present-day Hubble volume has, 
therefore, to be understood as averaging over the possible realization of the volume 
D inside the Hubble volume. This averaging is necessary if there is a big variance, 
which is however not expected, in the cosmic realizations of the volume D and has 
to be intended from a statistical point of view. 

However, this approach is not directly related to the averaged description the for- 
malism developed by Buchert is about. It should be seen as a possible interpretation 
of the latter. 
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Quantitative estimations 

Then, we have to check the strength and the properties of this effective pbr quan- 
titatively in order to see if it can explain the concordance model and how far it is 
from the standard cosmological constant. 

The following approaches have been mainly studied. A first possibility is by 
means of the perturbation theory. It has been found an instability in the perturba- 
tive expansion [36]. Since the perturbation approach breaks down, it is not possible 
to predict on firm grounds that backreaction is responsible for the present-day ac- 
celeration of the universe. However, it is intriguing that such an instability shows 
up only recently in the evolution of the universe. 

To overcome the problems coming from dealing perturbatively with non-linear 
perturbations, attention has been given to exact inhomogeneous solutions of Ein- 
stein's equations. 

In particular the Lemaitre-Tolman-Bondi (LTB) solution has been studied exten- 
sively in the literature [2l[8lia[3l[2lllini[6a[58l[62l[2l[35]. We will discuss this in 
Chapter [31 

Finally, numerical estimates were made in order to check quantitatively the im- 
portance of Qe). a plot of the absolute value of Qo, normalized by the global mean 
density, against scale for a Standard-CDM is shown in Fig. 12.51 The absolute values 
are shown there because the calculations were done in the newtonian approximation 
with periodic boundary conditions and therefore Qn is, by assumption, zero: Qd 
can indeed be written as a total divergence. The absolute value is then an estimate 
of the expected "backreaction" on some scale. 

From the plot we can see that the magnitude of the "backreaction" source term is of 
the same order as the mean density and higher on scales A < lOOMpc/h for SCDM. 
It quickly drops to a 10% effect on scales of A ~ 200Mpc/h. 

This is however only a hint toward a proper general relativistic estimation. The 
newtonian approximation neglects indeed important effects like, for example, the 
coupling fl2.26p between Qd and {TZ)^ without which there cannnot be a sizeable 
back-reaction. 
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Figure 2.5: The expected "backreaction" i.e., the quantity {QDl/^T^ip) d against the scale 

1/3 

L = a£) = Vjj (measured in Mpc/h) in linear /hnear (left panel) and log/log (right panel) 
format for a Standard-CDM model box of 1.8Gpc/h in size. This dimensionless quantity 
is still of the same order as the actual matter density on scales around lOOMpc/h jI3J /. 
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2.3 Smoothing out inhomogeneities - observational 
side 

The other possible approach focuses directly on the luminosity-distance-redshift 
relation d^lz) and looks at the effects inhomogeneities have on it. To carry out 
this analysis of the back-reaction it is necessary to have an inhomogeneous model 
to work with, in particular it is desirable to have an exact inhomogeneous solution 
of Einstein's equations. 

We remark that the dL{z) is a about the luminosity and redshift of photons that 
travelled through inhomogeneities and have their effects "integrated out" . In other 
words, even if we are not averaging explicitly the universe, the luminosity-distance- 
redshift relation has an implicit meaning of average in itself. 

This is the reason why this approach can be regarded as a possible answer to the 
smoothing proceses as well as the approches of the previous section are. 

It has been shown that the LTB solution can be used to fit the observed d^lz) 
without the need for dark energy (for example in [31 [Tj ) . To achieve this result , how- 
ever, it is necessary to place the observer at the center of a rather big underdensity. 
This is in contrast with the experimental verifications of the cosmological principle. 
To overcome this fine-tuning problem, in ChapterlUa Swiss-cheese model will be used 
where the observer is in the cheese and looks through Swiss-cheese holes constructed 
out of an LTB solution. 

It is of great interest to compare this physically meaningful approach with the 
light-cone average fitting we will use. Both are indeed tied to observations. 



Chapter 3 



The backreaction effect through a 
"homogeneous" dust LTB model 



In this chapter we will try to build an inhomogeneous universe model searching for 
a back-reaction effect of sub-Hubble perturbations. In particular we will aim at 
finding the average dynamics effects described in Sect. 12.2.11 

We will make use of the dust LTB model, described in Appendix [Bl It has been 
shown that the LTB solution can be used to fit the observed diiz) without the need 
of dark energy (for example in |3j). To achieve this result, however, it is necessary 
to place the observer at the center of a rather big underdensity. This is in contrast 
with the experimental verifications of the cosmological principle. 
To overcome this fine-tuning problem we will build a model that does not single out 
the center as a preferred place. To this end, we paid attention to choose appropriate 
densities and curvature profiles. 

To be clear, the density profile we have in mind is: 



In fig. l3.1l there is a sketch of that density profile. This shape is intermediate between 
a cubic lattice without preferred positions and a standard LTB model in which the 
center is singled out by its different density and/or spatial curvature. It is a sort 
of exact "spherical lattice": even though the inhomogeneities of the density are 
centered on the center of symmetry, an observer moving around such a universe 
should not note sizeable differences between one place and the other. 




(3.1) 
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Figure 3.1: An example of the density proGle l\3.1\) we are going to use in the Leniaitre- 
Tolman-Bondi modeL 



3.1 Flat case 

We will now consider a flat LTB model, 
with more elaborate models with E 0. 



This analysis will be useful when dealing 



3.1.1 Our model 

The density we will use has only inhomogeneities of one size: 

p(r, t) = Pm + Pm cos k r (3.2) 

and does not have the weak singularity mentioned in [64J. This density profile 
is clearly nothing more than a toy-model: the present-day universe is far more 
inhomogeneous. However, it will have the main features in order to look for a back- 
reaction effect. 

In f l3.2p k = 2t{/\ and we will choose A = 0.05 which, with the units of Table [BTTI 
corresponds to 420 Mpc. The voids will be roughly 20 times smaller than r^Bi the 
comoving distance traveled by a photon since the big bang. 

As said in Appendix [HI we fixed the gauge freedom by choosing i{r) = i and 
Y{r, t) = r. 



3.1. FLAT CASE 
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We will show how its density evolves in the next section when we will study the 
curved case. In the flat case the evolution is indeed reversed: structures are not 
forming, but spreading out. We will show better this in the Chapter HI 



3.1.2 Back-reaction effects 

We built this model in search of the back- reaction effects described in Sect. 12.2. It 
we are interested in the quantities on the right hand side of eq. f l2.2ip . which in the 
units of Table IB. II is: 

Hh = l{p)D~l{n)^-^Qn (3.3) 

We have given the analytical expression of their constituents for a LTB model in 
Appendix [Bl 

In the fiat case one finds out immediately from (IB. Ill) that (TZ) d = and with a 
short calculation that = 0, thus there is no bacicreaction at all i.e. the expansion 
rate of the LTB model and the one of the FRW model are the same. It exactly holds: 

QD = 1/2 (3.4) 

These results do not depend on the choice of the free functions of the LTB model. 
This is not trivial: in the fiat case the local equations are invariant under the 
averaging process: this is the cause of the absence of baclc-reaction. 

However is not possible to conclude that there is no baclcreaction in the curved 
case. In the fiat case we have, from eq. fl4.22p . = Y and this can be seen, 
alternatively, as the reason why we found no back-reaction: the equation for Y are 
already averaged as pointed out in Appendix [B] and therefore there can not be any 
different average dynamics. In the curved case, however, an Y. The expansion is 
now: 

If L^drOY^Y'/W Vd do , , 

VdJd j^drY^Y'/W Vd an 

where R marks the outer shell of the spherical domain Vd- It is not anymore true 
that Hd = Y/Y. 
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3.2 Curved case 

Introducing the curvature enriches the model with respect to the flat case. We 
need to introduce the curvature for two reasons: the first is to have a more natural 
evolution and the second is to have a non trivial domain dynamics. 

3.2.1 Our model 

Except for the non-zero curvature, we will use the same model studied in the fiat 
case. As explained in Appendix [HI the LTB curvature E is connected to the FRW 
curvature k hj k = —2E/r'^. (It is not normalized to unity.) If £^ ~ const, away 
from the observer there will be no clash with the CMB measurement of k = 0. 
However, the center will be singled out in contrast with our idea of "homogeneous" 
model. To avoid this, when we will introduce the curvature, we will demand to have 
-E = on the average. 

We could think to use an oscillatory curvature, similarly to the behavior of the 
density. However there is too much freedom in the choice of E: it is, together with 
p and i, a free function of the model. 

E, however, is fixed by the initial conditions: requesting E = forced the velocity 
of the shells to compensate the inhomogeneities to have a fiat space. Thus we can 
thinlc that a more natural situation will happen if such a compensation does not take 
place. To have a realistic evolution, we demand therefore that there are no initial 
peculiar velocities at time t, that is, to have an initial expansion H independent of 
r. From Eq. (1B.9P this implies: 



The graph of E{r) chosen in this way is shown in Fig. 14.41 as you can see we have 
exactly an oscillating behavior without the arbitrariness of requesting it. As seen 
from the figure, the curvature E{r) is small compared with unity. Indeed, in many 
formulae = (1 + 2^)^/2 ^ 1 + E appears, therefore one should compare E with 
1. In spite of its smallness, the curvature will play a crucial role to allow a realistic 
evolution of structures, as we will see in the next section. 



E{r) = ^H'^^^{i)r' 



1 M(r) 



(3.6) 



3.2. CURVED CASE 
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3.2.2 The dynamics 

In Fig. 13.31 we show the evolution of p{r,t) for three times: t = i = —0.8 (the Big 
Bang is at tsB = — 1), ^ = —0.4, and t = (corresponding to today). As you can see, 
overdense regions start contracting and become thin shells (mimicking structures), 
while underdense regions become larger (mimicking voids), and eventually occupy 
most of the volume. 

The densities are normalized with respect to the average density {p)d — PFRwif) 
where the latter is the density of the FRW model which is at the initial time 
PFRwit) = {p{r,t))E,- The equalities are not exact because the oscillations in the 
density smooth out not exactly within a finite volume: this is due to the growing of 
the surface area with the respect to the radius. 

Moreover, remember that r is only a label for the shell whose Euclidean position 
at time t is Y{r,t). In Fig. I3.3l we have normalized Y{r,t) using rpRw = '^{''^,t)/(^{t) 
where a{t) is again the scale factor of the FRW model just introduced. 
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Figure 3.3: Behavior of the density profiles p{r, t) with respect to rpRw = 

times t = t = —0.8 (blue), t = —0.4 (green) and t = to = (red). The values of pFRw{i) 

are 1, 2.8, and 25, for t = 0, —0.4, —0.8, respectively. 



3.2. CURVED CASE 
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3.2.3 Back-reaction effects 

We will now look again at the quantities on the right hand side of eq. fl2.2ip . in 
particular, Qd and {71) d- 

We have found that Qd is negligible while this is not the case for {71) d- The 
reason is simple: the only difference with respect to the flat model in which both 
are zero is the curvature. While E appears explicitly in {7Z)d, in Qd it has only an 
indirect effect coming from the fact that now aD y^Y. 

From the integrability condition (12.261) . we can therefore deduce that, being 
Qo and {7Z)d decoupled, the latter will evolve like a standard curvature term: 
{7Z)^ ~ l/ajj. In other words, we will have wbr — —1/3. We have verified this 
numerically: in Fig. IB. 101 we plotted the time evolution of {7V)d- 

<R>D 
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t 



Figure 3.4: Time evolution for {TI)d- It evolves as a FRW curvature term: (JZ) ~ l/a| 
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3.3 Discussion 

We have applied the averaging scheme of Sect. 12.2.11 to a toy model that features 
the growing of structures. We looked quantitatively at the importance of the back- 
reaction sources Qd and {71)d- While we found a (very small) contribution from 
the average curvature, no appreciable effects were coming from Qd. 

We trace this to the exact absence of back-reaction for the flat LTB model. The 
ultimate reason is that the equation automatically performs a "euclidean" average 
with respect to the local quantity Y: averaging something twice will not, of course, 
give sizeable effects. We believe that this is caused by spherical symmetry (Birkhoff 
theorem) . 

In the curved case, the curvature enters only indirectly in Q^i and therefore we still 
do not have interesting domain effects. 

From this study we learn that LTB models with observers at the center are not 
interesting to study the formalism developed by Buchert. This could be, however, 
due to a poor use of them more than to their deficiency with respect to back-reaction 
studies in general. Having already discussed the possibility that the questions put 
by Buchert's formalism are not the observationally-meaningful questions, we will 
try to make a better use of the LTB model. Having seen that spherical symmetry 
suppresses back-reaction effects, we will try to go beyond this limitation building a 
Swiss-cheese (toy) model for the universe in the next Chapter. 

We will conclude acknowledging that recently and independently the reference 
[8J developed a similar density profile, using, however, a different curvature: we 
have used an oscillating-around-zero curvature which does not single out the center 
of symmetry. 



Chapter 4 
Swiss Cheese 



In this chapter (see also [lH HSj) we explore a cosmo logical toy model in order to 
attempt to understand the role of large-scale non-linear cosmic inhomogeneities in 
the interpretation of observable data. The model is based on a Swiss-cheese model, 
where the cheese consists of the usual Friedmann-Robertson- Walker (FRW) solution 
and the holes are constructed out of a Lemaitre-Tolman-Bondi (LTB) solution. We 
are focusing on a Swiss-cheese model because, even if it is made of spherical sym- 
metric holes, it is not a spherical symmetric model as a whole. It is a first step to 
go beyond spherical symmetry which will turn out to be the main limitation of LTB 
solutions. 

This model will turn out to be well-suited to study the dual point of view sketched 
in Sect. I1.2.3[ See [9l [TTl [12] for other works in this direction. 

4.1 Our LTB model 

We are going to study a Swiss-cheese model where the cheese consists of the usual 
Friedmann-Robertson- Walker solution and the spherically symmetric holes are con- 
structed from a Lemaitre-Tolman-Bondi solution. The particular FRW solution we 
will choose is a matter-dominated, spatially-flat solution, i.e., the Einstein-de Sitter 
(EdS) model. 

In this section we will describe the LTB model parameters we have chosen. We refer 
to Appendix [B] for notation, units and an introduction to the LTB models. 

First of all, for simplicity we will choose t{r) = t; i.e., specifying the initial 
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conditions for each shell at the same moment of time. 



1 




Figure 4.1: The Apollonian Gasket. 

We will now choose p{r,i) and W{r) in order to match the fiat FRW model at 
the boundary of the hole: i.e., at the boundary of the hole p has to match the 
FRW density and W{r) has to go to unity. A physical picture is that, given a FRW 
sphere, all the matter in the inner region is pushed to the border of the sphere 
while the quantity of matter inside the sphere does not change. With the density 
chosen in this way, an observer outside the hole will not feel the presence of the hole 
as far as local physics is concerned (this does not apply to global quantities, such 
the luminosity-distance-redshift relation for example). So the cheese is evolving 
as an FRW universe while the holes evolve differently. In this way we can imagine 
putting in the cheese as many holes as we want, even with different sizes and density 
profiles, and still have an exact solution of the Einstein equations (as long as there 
is no superposition among the holes and the correct matching is achieved). The 
limiting picture of this procedure is the Apollonian Gasket of Fig. 14.11 where all the 
possible holes are placed, and therefore the model has the strange property that it 
is FRW nowhere, but it behaves as an FRW model on the average. This idea was 
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first proposed by Einstein and Straus 
To be specific, we choose p(r, t) to be 

pir,t) = pFRwit) ir>rh), 

where e = 0.0025, n = 0.42, a = Vh/lO, tm = 0.037, A = 50.59, and pFRw{t} = 25. 
In Fig. 14.21 we plot this chosen Gaussian density profile. The hole ends at r/j = 0.042 
which iq^ 350 Mpc and roughly 25 times smaller than tbb- Note that this is not a 
very big bubble. But it is an almost empty region: in the interior the matter density 
is roughly 10^ times smaller than in the cheese. Our model consists of a sequence 
of up to five holes and the observer is looking through them. The idea, however, is 
that the universe is completely filled with these holes, which form a sort of lattice 
as shown in Fig. 14.31 In this way an observer at rest with respect to a comoving 
cheese-FRW observer will see an isotropic CMB along the two directions of sight 
shown in Fig. 14.31 




Figure 4.2: The densities p{r,t) (solid curve) and p{r,t) (dashed curve). Here, i = —0.8 
(recall tsB = The hole ends at = 0.042. The matching to the FKW solution is 

achieved as one can see from the plot of p{r,t). 



^To get this number from Table IBTTI vou need to multiply by a{to) ~ 2.92. 
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Figure 4.3: Sketch of our Swiss-cheese model. An observer at rest with respect to a 
comoving cheese-FRW observer will see an isotropic CMB along the two directions of 
sight marked with dotted red lines. Three possible positions for an observer are shown. 

It is useful to consider the velocity of a shell relative to the FRW background. 
We define 

Avsh{r, t) = aiTBir, t) - apRwit), (4.2) 

where aLTBir,t) = Y{r,t)/r. To have a realistic evolution, we demand that there 
are no initial peculiar velocities at time t, that is, to have an initial expansion H 
independent of r: Avsh{r,t) = 0. From Eq. ( IB. 91) this implies 

E{r) = lHi,„it)r'-l.^. (4.3) 

The graph of E{r) chosen in this way is shown in Fig. 14. 4[ As seen from the 
figure, the curvature E{r) is small compared with unity. Indeed, in many formulae 

appears, therefore one should compare E with 1. In spite 



4.1. OUR LTB MODEL 



41 



of its smallness, the curvature will play a crucial role to allow a realistic evolution 
of structures, as we will see in the next section. 

Also in Fig. 14.41 we graph k{r) = —2E{r)/r'^, which is the generalization of the 
factor k in the usual FRW models. (It is not normalized to unity.) As one can see, 
k{r) is very nearly constant in the empty region inside the hole. This is another 
way to see the reason for our choice of the curvature function: we want to have in 
the center an empty bubble dominated by negative curvature. 

It is important to note that the dynamics of the hole is scale- in dependent: small 
holes will evolve in the same way as big holes. To show this, we just have to express 
Eq. (IB. 61) with respect to a generic variable r = r/g where g fixes the scale. If we 
change g, i.e., if we scale the density profile, we will find the same scaled shape for 
k{r) and the same time evolution. This property is again due to spherical symmetry 
which frees the inner shells from the influence of the outer ones: We can think of a 
shell as an infinitesimal FRW solution and its behavior is scale independent because 
it is a homogeneous and isotropic solution. 




Figure 4.4: Curvature E{r) and k{r) necessary for the initial conditions of no peculiar 
velocities. 
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4.2 The dynamics 

Now we explore the dynamics of this Swiss-cheese model. As we have said, the cheese 
evolves as in the standard FRW model. Of course, inside the holes the evolution is 
different. This will become clear from the plots given below. 

We will discuss two illustrative flat case where E[r) = 0, and a curved 

case where E{r) is given by Eq. fl4.3p . We are really interested only in the second 
case because the first will turn out to be unrealistic. But the flat case is useful to 
understand the djTiamics. 

4.2.1 The flat case 

In Fig. 14.51 we show the evolution of Y{r, t) for the fiat case, E{r) = 0. In the figure 
Y{r,t) is plotted for three times: t = i = —0.8 (recall tBB = —1), t = —0.4, and 
t = (corresponding to today). 

From Fig. 14.51 it is clear that outside the hole, i.e., for r > Vh, Y{r,t) evolves as a 
FRW solution, Y{r,t) oc r. However, deep inside the hole where it is almost empty, 
there is no time evolution to Y{r, t): it is Minkowski space. Indeed, thanks to spher- 
ical symmetry, the outer shells do not influence the interior. If we place additional 
matter inside the empty space, it will start expanding as an FRW universe, but at a 
lower rate because of the lower density. It is interesting to point out that a photon 
passing the empty region will undergo no redshift: again, it is just Minkowski space. 

This counterintuitive behavior (empty regions expanding slowly) is due to the 
fact that the spatial curvature vanishes. This corresponds to an unrealistic choice of 
initial peculiar velocities. To see this we plot the peculiar velocity that an observer 
following a shell r has with respect to an FRW observer passing through that same 
spatial point. The result is also shown in Fig. 14.51 where it is seen that matter is 
escaping from the high density regions. This causes the evolution to be reversed as 
one can see in Fig. 14.51 from the density profile at different times: structures are not 
forming, but spreading out. 

Remember that r is only a label for the shell whose Euclidean position at time 
t is Y{r,t). In the plots of the energy density we have normalized Y{r,t) using 
rpRW = Y{r,t)/a(t). 
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r 



Figure 4.5: Behavior ofY{r, t) with respect to r, the pecuUar velocities v{r, t) with respect 
to r, and the density profiles p{r,t) with respect to rpRw = Y{r,t)/a{t), for the fiat case 
at times t = t = —0.8, t = —0.4 and t = Iq = 0. The straight lines for Y{r,t) are the 
FRW solutions while the dashed lines are the LTB solutions. For the peculiar velocities, 
matter is escaping from high density regions. The center has no peculiar velocity because 
of spherical symmetry, and the maximum of negative peculiar velocity is before the peak 
in density. Finally, the values of p{oo,t) are 1, 2.8, and 25, for t = 0, —0.4, —0.8, 
respectively. 
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4.2.2 The curved case 

Now we will move to a more interesting and relevant case. We are going to use 
the E{r) given by Eq. (14. 3p : the other parameters will stay the same. Comparison 
with the flat case is useful to understand how the model behaves, and in particular 
the role of the curvature. In Fig. 14.61 the results for Y{r,t) in the curved case are 
plotted. Time goes from t = t = —0.8 to t = 0. 

As one can see, now the inner almost empty region is expanding faster than the 
outer (cheese) region. This is shown clearly in Fig. 14.71 where also the evolution of 
the inner and outer sizes is shown. Now the density ratio between the cheese and the 
interior region of the hole increases by a factor of 2 between t = i and t = 0. Initially 
the density ratio was 10^, but the model is not sensitive to this number since the 
evolution in the interior region is dominated by the curvature (fc(r) is much larger 
than the matter density). We stress now the fact that the crucial ingredient is to 
have a faster-than-cheese expanding void. 

The peculiar velocities are now natural: as can be seen from Fig. 14.61 matter 
is falling towards the peak in the density. The evolution is now realistic, as one 
can see from Fig. 14. 6[ which shows the density profile at different times. Overdense 
regions start contracting and they become thin shells (mimicking structures), while 
underdense regions become larger (mimicking voids), and eventually they occupy 
most of the volume. 

Let us explain why the high density shell forms and the nature of the shell cross- 
ing. Because of the distribution of matter, the inner part of the hole is expanding 
faster than the cheese; between these two regions there is the initial overdensity. It 
is because of this that there is less matter in the interior part. (Remember that we 
matched the FRW density at the end of the hole.) Now we clearly see what is hap- 
pening: the overdense region is squeezed by the interior and exterior regions which 
act as a clamp. Shell crossing eventually happens when more shells - each labeled 
by its own r - are so squeezed that they occupy the same position Y, i.e. when 
y = 0. Nothing happens to the photons other than passing through more shells at 
the same time: this is the meaning of the Qrr metric coefficient going to zero. 

A remark is in order here: In the inner part of the hole there is almost no matter, 
it is empty. Therefore it has only negative curvature, which is largely dominant over 
the matter: it is close to a Milne universe. 
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CURVED CASE 




r 



Figure 4.6: Behavior ofY{r, t) with respect to r, the pecuUar velocities v{r, t) with respect 
to r, and the density profiles p{r,t) with respect to rpRw = Y{r,t)/a{t), for the curved 
case at times t = t = —0.8, t = —0.4 and t = to = 0. The straight lines for Y{r, t) are the 
FEW solutions while the dashed lines are the LTB solutions. For the peculiar velocities, 
the matter gradually starts to move toward high density regions. The solid vertical line 
marks the position of the peak in the density with respect to r. For the densities, note 
that the curve for p{r,0) has been divided by 10. Finally, the values of p{oo,t) are 1, 2.8, 
and 25, for t = 0, —0.4, —0.8, respectively. 
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Figure 4.7: Evolution of the expansion rate and the size for the inner and outer regions. 
Here "inner" refers to a point deep inside the hole, and "outer" refers to a point in the 
cheese. 
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4.3 Photons 

We are mostly interested in observables associated with the propagation of photons 
in our Swiss-cheese model: indeed, our aim is to calculate the luminosity-distance- 
redshift relation dL{z) in order to understand the effects of inhomogeneities on ob- 
servables. Our setup is illustrated in Fig. 14.81 where there is a sketch of the model 
with only 3 holes for the sake of clarity. Notice that photons are propagating through 
the centers. 

We will discuss two categories of cases: 1) when the observer is just outside the 
last hole as in Fig. 14.81 and 2) when the observer is inside the hole. The observer 
in the hole will have two sub ) the observer located on a high-density shell, 

and b) the observer in the center of the hole. We are mostly interested in the first 
case: the observer is still a usual FRW observer, but looking through the holes in 
the Swiss cheese. 




-Fh Th 3rh 5rh 



Figure 4.8: Sketch of our model in comoving coordinates. The shading represents the 
initial density profile: darker shading implies larger densities. The uniform gray is the 
FRW cheese. The photons pass through the holes as shown by the arrow. 
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4.3.1 Finding the photon path: an observer in the cheese 

We will discuss now the equations we will use to find the path of a photon through 
the Swiss cheese. The geodesic equations can be reduced to a set of four first-order 
differential equations (we are in the plane 6 = 7r/2): 



— = — -I (z + ir-^ \ -d^ 



dz 
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where A is an affine parameter that grows with time. The third equation is actually 
the null condition for the geodesic. Thanks to the initial conditions chosen we have 
,2(0) = 0. These equations describe the general path of a photon. To solve the 
equations we need to specify the constant c^, a sort of angular momentum density. 
A first observation is that setting c,^ = allows us to recover the equations that 
describe a photon passing radially trough the centers: dt/dr = Y'/W. 

We are interested in photons that hit the observer at an angle a and are passing 
trough all the holes as shown in Fig. 14.81 To do this we must compute the inner 
product of X* and y^, which are the normalized spatial vectors tangent to the radial 
axis and the geodesic as shown in Fig. 14.91 A similar approach was used in Ref . [2] . 
The inner product of and is expressed through 
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(4.7) 
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The vectors are anchored to the shell labeled by the value of the affine parameter 
A = 0, that is, to the border of the hole. Therefore, they are relative to the comoving 
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Figure 4.9: A photon hitting the observer at an angle a. 



observer located there. In the second equation we have used the initial conditions 
given in the previous set of equations, while to find the last equation we have used 
the null condition evaluated at A = 0. 

The above calculations use coordinates relative to the center. However, the angle 
a is a scalar in the hypersurface we have chosen: we are using the synchronous and 
comoving gauge. Therefore, a is the same angle measured by a comoving observer 
of Fig. 14.91 located on the shell r = — r^: it is a coordinate transformation within 
the same hypersurface. 

Given an angle a we can solve the equations. We have to change the sign in Eq. 
(14.41) when the photon is approaching the center with respect to the previous case 
where it is moving away. Also, we have to sew together the solutions between one 
hole and another, giving not only the right initial conditions, but also the appropriate 
constants c<^ (see Appendix \C\i . 

Eventually we end up with the solution t(A), t(A), 0(A) and z{X) from which we 
can calculate the observables of interest. 
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4.3.2 Finding the photon path: an observer in the hole 

Finding the solution in this case is the same as in the previous case with the only 
difference that in Eq. fl4.4p the initial condition is now r(0) = Vobs- But this observer 
has a peculiar velocity with respect to an FRW observer passing by. This, for 
example, will make the observer see an anisotropic cosmic microwave background 
as it is clear from Fig. I4.3[ This Doppler effect, however, is already corrected in the 
solution we are going to find since we have chosen z{0) = as initial condition. 

There is however also the effect of light aberration which changes the angle a seen 
by the comoving observer with respect to the angle ap^w seen by an FRW observer. 
The photon can be thought as coming from a source very close to the comoving 
observer: therefore there is no peculiar motion between them. The FRW observer 
is instead moving with respect to this reference frame as pictured in Fig. 14.101 The 
relation between a and apuw is given by the relativistic aberration formula: 

cosa + v/c 

cos apnw = — — -, • (4.9) 

1 + v/c cos a 

The angle changes because the hypersurface has been changed. The velocity will be 
taken from the calculation (see Fig. 14.61 for the magnitude of the effect). 




Figure 4.10: A comoving observer and a FRW observer live in different frames, this results 
in a relative velocity vpRW between observers. 
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4.3.3 Distances 



The angular diameter distance is defined as: 



dA 



D 



(4.10) 



C^FRW 



where D is the proper diameter of the source and a is the angle at which the source 
is seen by the observer. Using this definition to find (Ia we have 



The formula we are going to use for dA is exact in the limit of zero curvature. 
However in our model E{r) is on average less than 0.3% and never more than 0.4%, 
as it can be seen from Fig. 14. 4t therefore the approximation is good. Moreover, we 
are interested mainly in the case when the source is out of the last hole as pictured 
in Fig. 14.81 and in this case the curvature is exactly zero and the result is exact. 

We have checked that the computation of dA is independent of a for small an- 
gles and that the result using the usual FRW equation coincides with theoretical 
prediction for dA- We also checked that dA reduces to Y{r,t) when the observer is 
in the center. 

Finally we checked our procedure in comparison with the formula {E.31) of Ref. 
[8]: this is a rather different way to find the angular distance and therefore this 
agreement serves as a consistency check. We placed the observer in the same way 
and we found the same results provided that we use the angle a uncorrected for the 
light-aberration effect. 

4.4 Results: observer in the cheese 

Now we will look through the Swiss cheese comparing the results with respect to a 
FRW-EdS universe and a ACDM case. 

We will first analyze in detail the model with five holes, which is the one which 
we are most interested in. For comparison, we will study models with one big hole 



dA 




(4.11) 



The luminosity distance will then be: 



dL={l + zfdA. 



(4.12) 
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and one small hole. In the model with one big hole, the hole will be five-times bigger 
in size than in the model with five holes: i.e., they will cover the same piece of the 
universe. 

The observables on which we will focus are the changes in redshift z{X), angular- 
diameter distance dA^z), luminosity distance d^iz), and the corresponding distance 
modulus A'm{z). 

4.4.1 Redshift histories 

Now we will first compare the redshift undergone by photons that travel through 
the model with either five holes or one hole to the FRW solution of the cheese. In 
Fig. 14. Ill the results are shown for a photon passing through the center with respect 
to the coordinate radius. As one can see, the effects of the inhomogeneities on the 
redshift are smaller in the five-hole case. 

It is natural to expect a compensation, due to the spherical symmetry, between 
the ingoing path and the outgoing one inside the same hole. This compensation is 
evident in Fig. 14. Ill 

However, there is a compensation already on the scale of half a hole as it is clear 
from the plots. This mechanism is due to the density profile chosen, that is one 
whose average matches the FRW density of the cheese: roughly speaking we know 
that z' = H oc p = Pphw + ^P- We chose the density profile in order to have {5p) = 0, 
and therefore in its journey from the center to the border of the hole the photon will 
see a (H) ~ Hpi^w and therefore there will be compensation for z'. It is somewhat 
similar to the screening among positive and negative charges. 

Let us see this analytically. We are interested in computing a line average of the 
expansion along the photon path in order to track what is going on. Therefore, we 
shall not use the complete expansion scalar: 

= T'^, = 2y + ^, (4.13) 
but, instead, only the part of it pertinent to a radial line average: 

e,, = Tl, = ^^Hr, (4.14) 
where Fq^, are the Christoffel symbols and 6 is the trace of the extrinsic curvature. 
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Figure 4.11: Redshift histories for a photon that travels from one side of the one-hole 
chain (left) and five-hole chain (right) to the other where the observer will detect it at 
present time. The "regular" curve is for the FRW model. The vertical lines mark the 
edges of the holes. The plots are with respect to the coordinate radius r. Notice also that 
along the voids the redshift is increasing faster: indeed z'(r) = H{z) and the voids are 
expanding faster. 



Using Hr, we obtain: 

, ^ f"" dr Hr Y' /W Y 

where the approximation comes from neglecting the (small) curvature and the last 
equality holds thanks to the density profile chosen. This is exactly the result we 
wanted to find. However, we have performed an average at constant time and there- 
fore we did not let the hole and its structures evolve while the photon is passing: 
this effect will partially break the compensation. This sheds light on the fact that 
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photon physics seems to be affected by the evolution of inhomogeneities more than 
by inhomogeneities themselves. We can argue that there should be perfect compen- 
sation if the hole will have a static metric such as the Schwarzschild one. In the end, 
this is a limitation of our assumption of spherical symmetry. 

This compensation is almost perfect in the five-hole case, while it is not in the 
one-hole case: in the latter case the evolution has more time to change the hole while 
the photon is passing. Summarizing, the compensation is working on the scale Vh of 
half a hole. These results are in agreement with Ref. [S]. 

From the plot of the redshift one can see that function z{r) is not monotonic. 
This happens at recent times when the high-density thin shell forms. This blueshift 
is due to the peculiar movement of the matter that is forming the shell. This feature 
is shown in Fig. 14.121 where the distance between the observer located just out of 
the hole at r = and two different shells is plotted. In the solid curve one can see 
the behavior with respect to a normal redshifted shell, while in the dashed curve 
one can see the behavior with respect to a shell that will be blueshifted: initially 
the distance increases following the Hubble flow, but when the shell starts forming, 
the peculiar motion prevails on the Hubble flow and the distance decreases during 
the collapse. 

It is finally interesting to interpret the redshift that a photon undergoes passing 
the inner void. The small amount of matter is subdominant with respect to the 
curvature which is governing the evolution, but still it is important to define the 
space: in the limit of zero matter in the interior of the hole, we recover a Milne 
universe, which is just (half of) Minkowski space in unusual coordinates. Before 
this limit the redshift was conceptually due to the expansion of the spacetime, after 
this limit it is instead due to the peculiar motion of the shells which now carry no 
matter: it is a Doppler effect. 

4.4.2 Luminosity and Angular-Diameter Distances 
The five-hole model 

In Fig. 14.131 the results for the luminosity distance and angular distance are shown. 
The solution is compared to the one of the ACDM model with Qm = 0.6 and 
^DE = 0.4. Therefore, we have an effective go = ^a//2 — ^Ide = —0.1. In all the 
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Figure 4.12: Distance between the observer and two different sliells. In the solid curve 
r = 0.55 r/j will be redshifted, while in the dashed curve, r = 0.8 r/j will be blueshifted. 
The latter indeed will start to collapse toward the observer. Time goes from t = —0.8 to 
t = 0. The observer is located just outside of the hole at r = r^. 



plots we will compare this ACDM solution to our Swiss-cheese solution. The strange 
features which appear near the contact region of the holes at recent times are due 
to the non-monotonic behavior of z{r), which was explained in the previous section. 

The distance modulus is plotted in the top panel of Fig. 14. 131 The solution shows 
an oscillating behavior which is due to the simplification of this toy model in which 
all the voids are concentrated inside the holes and all the structures are in thin 
spherical shells. For this reason a fitting curve was plotted: it is passing through the 
points of the photon path that are in the cheese between the holes. Indeed, they are 
points of average behavior and well represent the coarse graining of this oscillating 
curve. The simplification of this model also tells us that the most interesting part 
of the plot is farthest from the observer, let us say aX z > 1. In this region we can 
see the effect of the holes clearly: they move the curve from the EdS solution (in 
purple) to the ACDM one with = 0.6 and Qde = 0.4 (in blue). Of course, the 
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Figure 4.13: On the bottom the luminosity distance dL{z) in the Rve-hole model (jagged 
curve) and the ACDM solution with Qm = 0.6 and ^Ide = 0.4 (regular curve) are shown. 
In the middle is the change in the angular diameter distance, A(i^(z), compared to a 
ACDM model with Qm = 0.6 and VLqe = 0.4. The top panel shows the distance modulus 
in various cosmological models. The jagged line is for the Rve-hole LTB model. The 
regular curves, from top to bottom, are a ACDM model with i^M = 0.3 and ^de = 0.7, a 
ACDM model with Qm = 0.6 and ^^de = 0.4, the best smooth fit to the LTB model, and 
the EdS model. The vertical lines mark the edges of the five holes. 
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model in not realistic enough to reach the "concordance" solution. 

Here we are discussing a comparison of our results with those of Ref. In that 
paper they do not find the big difference from FRW results that we do. First of 
all, we can notice that we are able to reproduce their results using our techniques. 
The difference between their results and ours is that our model has very strong 
nonlinear evolution, in particular, close to shell crossing where we have to stop 
our calculations. Also, the authors of Ref. p] used smaller holes with a different 
density/initial- velocity profile. This demonstrated that a big change in observables 
may require either non-spherical inhomogeneities, or evolution very close to shell 
crossing. (We remind the reader that caustics are certainly expected to form in cold 
dark matter models.) 

Let us return now to the reason for our results. As we have seen previously, due 
to spherical symmetry there are no significant redshift effects in the five-hole case. 
Therefore, these effects must be due to changes in the angular-diameter distance. 
Fig. 14.141 is useful to understand what is going on: the angle from the observer is 
plotted. Through the inner void and the cheese the photon is going straight: they 
are both FRW solutions even if with different parameters. This is shown in the plot 
by constancy of the slope. The bending occurs near the peak in the density where 
the Qrr coefficient of the metric goes toward zero. Indeed the coordinate velocity of 
the photon can be split into an angular part: = d(f)/dt = 1 / ^ and a radial 
part Vr = dr / dt = \ j ^g^r- While behaves well near the peak, goes to infinity 
in the limit where shell crossing is reached: the photons are passing more and more 
matter shells in a short interval of time as the evolution approaches the shell-crossing 
point. Although in our model we do not reach shell crossing, this is the reason for 
the bending. We can, therefore, see that all the effects in this model, redshift and 
angular effects, are due to the evolution of inhomogeneities and this is primary due 
to the presence of a faster-than-cheese expanding void which, we think, is a crucial 
ingredient. 

The one-hole model: the big hole case 

Let us see now how the results change if instead of the five-hole model we use the 
one-hole model. We have already shown the redshift results in the previous section. 
As one can see from Fig. 14.151 the results are more dramatic: for high redshifts 
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Figure 4.14: The angle from the observer is plotted. The dashed vertical lines near the 
empty region mark the shell of maximum peculiar velocities of Fig. 14.61 The shaded regions 
represent the inner FRW solution. The solid vertical lines mark the peak in density. The 
angle at which the photon hits the observer is 2.7° on the left. 

the Swiss-cheese curve can be fit by a ACDM model with less dark energy than 
^DE = 0.6 as in the five-hole model. Nonetheless, the results have not changed 
so much compared to the change in the redshift effects discussed in the previous 
section. Indeed the compensation scale for angular effects is 2rh while the one for 
redshift effects is r^- 

The one-hole model: the small hole case 

Finally if we remove four holes from the five-hole model, we lose almost all the 
effects. This is shown in Fig. 14.161 now the model can be compared to a ACDM 
model with Qm = 0.95 and Qde = 0.05. 
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Figure 4.15: On the bottom is shown the luminosity distance diiz) in the one-hole model 
(jagged curve) and the ACDM solution with CIm = 0.6 and ^de = 0.4 (regular curve). 
In the middle is the change in the the angular diameter distance, AdA{z), compared to a 
ACDM model with = 0.6 and Q,£,e = 0.4. On the top is shown the distance modulus 
in various cosmological models. The jagged line is for the one-hole LTB model. The 
regular curves, from top to bottom are a ACDM model with f^Af = 0.3 and Qde = 0.7, a 
ACDM model with Qm = 0.6 and Qde = 0.4 and the EdS model. The vertical lines mark 
the edges of the hole. 
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Figure 4.16: On the bottom is shown the luminosity distance di^z) in the 1-hole model 
(jagged curve) and the ACDM solution with Qm = 0.95 and CIde = 0.05 (regular curve). 
In the middle is the change in the the angular diameter distance, AdA{z), compared to 
a ACDM model with Q^j = 0.95 and VLde = 0.05. On the top is shown the distance 
modulus in various cosmological models. The jagged line is for the one-hole LTB model. 
The regular curves, from top to bottom are a ACDM model with Um = 0.3 and Qde = 0.7, 
a ACDM model with Qm = 0.95 and ^de = 0.05 and the EdS model. The vertical lines 
mark the edges of the hole. 
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4.5 Results: observer in the hole 

Now we will examine the case in which the observer is inside the last hole in the 
five-hole model. We will first put the observer on the high-density shell and then 
place the observer in the center. 

4.5.1 Observer on the high density shell 

In the section we show the results for the observer on the high-density shell. As 
one can see from Fig. 14.171 now the compensation in the redshift effect is lost: the 
photon is not completing the entire last half of the last hole. The results for the 
luminosity distance and the angular distance do not change much as shown in Fig. 

Km 

Remember that in this case the observer has a peculiar velocity compared to the 
FRW observer passing through the same point. We correct the results taking into 
account both the Doppler effect and the light aberration effect. 




Figure 4.17: Redshift histories for a photon that travels through the five-hole-chain to 
the observer placed on the high density shell. The "regular" line is for the FRW model. 
A is the afRne parameter and it grows with the time which go from the left to the right. 
The vertical lines mark the end and the beginning of the holes. 
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Figure 4.18: On the bottom is shown the luminosity distance diiz) in the Eve-hole model 
(jagged curve) and the ACDM solution with CIm = 0.6 and ^Ide = 0.4 (regular curve). In 
the middle is the change in the angular diameter distance, AdA(z), compared to a ACDM 
model with Qm = 0.6 and ^de = 0.4. On the top is shown the distance modulus in 
various cosmological models. The jagged line is for the Eve-hole LTB model. The regular 
curves, from top to bottom are a ACDM model with Qm = 0.3 and i^^E = 0.7, a ACDM 
model with = 0.6 and Qde = 0.4, the best smooth Et to the LTB model, and the EdS 
model. The vertical lines mark the edges of the Eve holes. 
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4.5.2 Observer in the center 

In this section we show the results for the observer in the center. As confirmed 
by Fig. 14.191 the compensation in the redshift effect is good: the photon is passing 
through an integer number of half holes. 

The results for the luminosity distance and the angular distance look worse as 
shown in Fig. 14.201 but this is mainly due to the fact that now the photon crosses 
half a hole less than in the previous cases and therefore it undergoes less bending. 

In this case the observer has no peculiar velocity compared to the FRW one: this 
is a result of spherical symmetry. 




Figure 4.19: Redshift histories for a photon that travels through the five-hole-chain to 
the observer placed in the center. The "regular" line is for the FRW model. A is the affine 
parameter and it grows with the time which go from the left to the right. The vertical 
lines mark the end and the beginning of the holes. 
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Figure 4.20: The bottom panel shows the luminosity distance dL{z) in the Rve-hole model 
(jagged curve) and the ACDM solution with CIm = 0.6 and ^Ide = 0.4 (regular curve). In 
the middle is the change in the angular diameter distance, AdA(z), compared to a ACDM 
model with CIm = 0.6 and Qde = 0.4. On the top panel the distance modulus in various 
cosmological models is shown. The jagged line is for the Eve-hole LTB model. The regular 
curves, from top to bottom are a ACDM model with Qm = 0.3 and ^^de = 0.7, a ACDM 
model with = 0.6 and Qde = 0.4, the best smooth fit to the LTB model, and the EdS 
model. The vertical lines mark the edges of the five holes. 
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4.6 The fitting problem 

Now that we have seen how the luminosity-distance-redshift relation is affected by 
inhomogeneities, we want to study the same model from the point of view of light- 
cone averaging to see if we can gain insights into how inhomogeneities renormalize 
the matter Swiss-cheese model and mimic a dark-energy component. 

We are going to work out in this section the ideas introduced in Section I2.2.2[ 
A remark is, however, in order here: in the previous section we did not fit the diiz) 
with an FRW solution. We have simply compared the shape of the diiz) for the 
Swiss-cheese model with the one of a ACDM model. 

We intend now to fit a phenomenological FRW model to our Swiss-cheese model. 
The FRW model we have in mind is a spatially fiat model with a matter component 
with present fraction of the energy density Qm = 0.25, and with a phenomenological 
dark-energy component with present fraction of the energy density Q\ = 0.75. We 
will assume that the dark-energy component has an equation of state 

w{a) =Wo + Wa ll - — ] = Wo + Wa (4.16) 
\ aoj l + z 

Thus, the total energy density in the phenomenological model evolves as 

^ = Qm{1 + Zf + fiA(l + ^)3(l+-o+-.) (^^^^\ . (4.17) 

We will refer to this model as the phenomenological model. 

Our Swiss-cheese model is a lattice of holes as sketched in Fig. 14. 8[ the scale of 
inhomogeneities is therefore simply the size of a hole. We are interested in under- 
standing how the equation of state of "dark energy" in the phenomenological model 
changes with respect to r^, and in particular, why. Of course, in the limit r/^ — 0, 
we expect to find w = 0, that is, the underlying EdS model out of which the cheese 
is constructed. 
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The procedure developed by Ref. [22] is summarized by Fig. I4.21[ We refer the 
reader to that reference for a more thorough analysis and to Ref. [22] and references 
therein for recent developments. We will focus now in using our Swiss-cheese model 
as a cosmological (toy) model. 




Figure 4.21: In the null data best fitting, one successively chooses maps from the real 
cosmological model U to the FRW model U' of the null cone vertex p' , the matter 4- 
velocity at p' , a two-sphere S' on the null cone of p' and a point q' on the 2-sphere. This 
establishes the correspondence ■0 of points on the past null cone ofp', C^{p'), to the past 
null cone of p, C~{p'), and then compares initial data at q' and at q. From Figure 2 of 



4.6.1 Choice of vertex points 

We start choosing the two observers to be compared. In the homogeneous FRW 
model every observer is the same thanks to spatial homogeneity. We choose an 
observer in the cheese as the corresponding observer in our Swiss-cheese model, in 
particular the one shown in Fig. 14. 8[ 

Our model allows us to choose also the time of observation, which, in general, is 
a final product of the comparison. We now explain why. 

The FRW model we will obtain from the fit will evolve differently from the 
Swiss cheese: the latter evolves as an EdS model, while the former will evolve as 
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a quintessence-like model. They are really different models. They will agree only 
along the light cone, that is, on our observations. 

Now, for consistency, when we make local measurementj^ the two models have 
to give us the same answer: local measurements indeed can be seen as averaging 
measurements with a small enough scale of averaging, and the two models agree 
along the past light cone. 

Therefore, we choose the time in order that the two observers measure the same 
local density. This feature is already inherent in Eq. (14.1 7p : the phenomenological 
model and the Swiss-cheese model evolve in order to have the same local density, 
and therefore the same Hubble parameter, at the present time. 

4.6.2 Fitting the 4-velocity 

The next step is to fit the four-velocities of the observers. In the FRW model we 
will choose a comoving observer, the only one who experiences an isotropic CMB. 
In the Swiss-cheese model, we will choose, for the same reason, a cheese-comoving 
observer. Again, our Swiss-cheese model considerably simplifies our work. 

4.6.3 Choice of comparison points on the null cones 

Now that the past null cones are uniquely determined, we have to choose a measure 
of distance to compare points along each null cone. 

First, let us point out that instead of the entire two-sphere along the null cone, we 
will examine, only a point on it. This is because of the simplified set-up of our Swiss- 
cheese model in which the observer is observing only in two opposite directions, as 
illustrated in Fig. 14. 8[ This means that we can skip the step consisting in averaging 
our observable quantities over the surface of constant redshift, which is generally 
necessary in order to be able to compare an inhomogeneous model with the FRW 
model PJ. 

^Conceptually, it could not be possible with a realistic universe model to make local mea- 
surements that could be directly compared to the smoothed FRW model. We are allowed to do 
so thanks to our particular Swiss-cheese model in which the cheese well represents the average 
properties of the model. 
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Coming back to the main issue of this section, we will use the observed redshift 
z to compare points along the null cones. Generally, the disadvantage of using it 
is that it does not directly represent distances along the null cone. Rather, the 
observed value z is related to the cosmological redshift zc by the relation: 

l + z^{l + zo){l + zc){l + zs) (4.18) 

where zq is the redshift due to the peculiar velocity of the observer O and Zs that 
due to the peculiar velocity of the source. The latter, in particular, is a problem 
because local observations cannot distinguish Zg from Zc- 

However, our set up again simplifies this task. The chosen observers are, indeed, 
both comoving (in the Swiss-cheese model because the observer is in the cheese, and 
in the phenomenological model by construction), and therefore zq = 0. Regarding 
the sources, we know exactly their behavior because we have a model to work with. 

The sources are also comoving; however, there are structure-formation effects 
that should be disentangled from the average evolution. For this reason we will 
perform averages between points in the cheese (the meaning of this will be clear in 
the next section) in order to smooth out these structure-formation effects. 
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Figure 4.22: An illustration of the points chosen for the averaging procedure. 



Now we are ready to set up the fitting of our Swiss-cfieese model. Ref. [31] stud- 
ied tfie approacli based on volume averaging outlined in Ref. [32] • This approach, 
however, is appropriate for studies concerning global dynamics, as in Refs. [TU| [T3] . 
As stressed previously, here we are instead interested in averages directly related to 
observational quantities, and we constructed our model following this idea: it is a 
model that is exactly solvable and "realistic" (even if still a toy) at the price of no 
interesting volume-averaged dynamics. 

Therefore, we will follow a slightly different approach from the ones outlined in 
Ref. [32j: we are going to fit averages along the light cone. This method will be 
intermediate between the fitting approach and the averaging approach. 

We will focus on the expansion scalar and the density. We will see that these two 
quantities behave differently under averaging. We denote by Q^'"{r,t) a quantity in 
the Swiss-cheese model we want to average. We denote by Q^^'^it) the corresponding 
quantity we want to fit to the average of Q^'^{r,t). Note that Q^^^(t) does not 
depend on r because the phenomenological model we will employ to describe the 
Swiss-cheese model is homogeneous. 

Again, the fit model is a phenomenological homogeneous model (just refereed to 
as the phenomenological model). It need not be the model of the cheese. 

The procedure is as follows. First we will average (5^^(?", t) for a photon that 
starts from the emission point E of the five-hole chain and arrives at the locations of 
observers Oi of Fig. 14.221 We have chosen those points because they well represent 
the average dynamics of the model. Indeed, these points are not affected by structure 
evolution because they are in the cheese. Then, we will compare this result with the 
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average of Q^^'^{t) for the phenomenological and homogeneous source with density 
given by Eq. (14.1 7p with an equation of state w given by Eq. (14. 161) . 
The two quantities to be compared are therefore: 



(g'"^)^. = / drY'/W 
'e 



dr Q^^{r,t{r)) Y'{r,t{r))/W{r) 

E 



dr a. 



FIT 
E 



1 rO, 



CbpjT (W(r)), (4.19) 



where t{r) and tpir^r) are the photon geodesies in the Swiss-cheese model and in the 
phenomenological one, respectively. The functions tprr^r), apjT and other quantities 
we will need are obtained solving the Friedman equations with a source described 
by Eq. (14.171) with no curvature. The points Oi in the Swiss-cheese model of Fig. 
14.221 are associated to points in the phenomenological model with the same redshift, 
as discussed in Sec. 14.6.31 

We will then find the w that gives the best fit between {Q^^'^) and (Q^^), that 
is, the choice that minimizes: 

(4-20) 

i 

Of course, in the absence of inhomogeneities, this method would give w = 0. 
Let us summarize the approach: 

• We choose a phenomenological quintessence-like model that, at the present 
time, has the same density and Hubble parameter as the EdS-cheese model. 

• We make this phenomenological model and the Swiss-cheese model correspond 
along the light cone via light-cone averages of Q. 

• We can substitute the Swiss-cheese model with the phenomenological model 
as far as the averaged quantity Q is concerned. 

The ultimate question is if it is observationally meaningful to consider Q, as op- 
posed to the other choice of domain averaging at constant time, which is not directly 
related to observations. We will come back to this issue after having obtained the 
results. 
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Averaged expansion 

The first quantity in wliich we are interested is the expansion rate. To average 
the expansion rate we will follow the formalism developed in Sec. I4.4.1I We will 
therefore apply Eqs. (14.19114 .201) to Q^'^ = Hr = Y' /Y', where we remember that Hr 
is the radial expansion rate. The corresponding quantity in the phenomenological 
model is Q^^^ = apu-fapiT- 

For the same reason there is good compensation in redshift effects (see Sec. 
14.4.11) . we expect (Hr) to behave very similarly to the FRW cheese solution. Indeed, 
as one can see in Fig. 14.231 the best fit of the Swiss-cheese model is given by a 
phenomenological source with w ~ 0, that is, the phenomenological model is the 
cheese-FRW solution itself as far as the expansion rate is concerned. 




Figure 4.23: Average expansion rate. The yellow points are {H^^)^. while the crosses 
are {H^^'^)^. . EOi means that the average was performed from E and Oi with respect 
to Fig. 14.221 The best fit is found for w ~ 0, that is, the phenomenological model is the 
cheese-FRW solution itself as far as the expansion rate is concerned. 
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Averaged density 

The situation for tlie density is very different. The photon is spending more and 
more time in the (large) voids than in the (thin) high density structures. We apply 
Eqs. f l4.19lH720l) to Q^*" = p^*". The corresponding quantity in the phenomenological 
model is Q^^^ = p^^^ where p^^""" is given by Eq. fl4.17p . The results are illustrated 
in Fig. I4.24t the best fit is for wq = -1.95 and Wa = 4.28. 




Figure 4.24: Average density in pco units. The yellow points are {p while the 

crosses are {p^^'^)^^- EOi means that the average was performed from E and Oi with 
respect to Fig. \4.22[ The parametrization of p^^^ is from Eq. i4.17\) . The best fit is found 
for wq = —1.95 and Wa = 4.28. 

As we will see in Sec. 14. 7[ we can achieve a better fit to the concordance model 
with smaller holes than the ones of 350 Mpc considered here. We anticipate that 
for a holes of radius = 250 Mpc, we have wq = —1.03 and Wa = 2.19. 

We see, therefore, that this Swiss-cheese model could be interpreted, in the 
FRW hypothesis, as a homogeneous model that is initially dominated by matter 
and subsequently by dark energy: this is what the concordance model suggests. We 
stress that this holds only for the light-cone averages of the density. 
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4.6.5 Explanation 

Let us first explore the basis for what we found. In Fig. 14.251 we show the density 
along the light cone for both the Swiss-cheese model and the EdS model for the 
cheese. It is clear that the photon is spending more and more time in the (large) 
voids than in the (thin) high density structures. 




Figure 4.25: Density along the light-cone for the Swiss-cheese model (the spiky curve) 
and the EdS model of the cheese (the regular curve). The labeling of the x-axis is the 
same one of Fig. \4.22[ 



To better show this, we plotted in Fig. 14.261 the constant-time, line-averaged 
density as a function of time. The formula used for the Swiss-cheese model is 

dr p{r,t)Y'{r,t)/W{r) / j drY'/W , (4.21) 

while for the cheese, because of homogeneity we can just use p(t) of the EdS model. 
As one can see, the photon is encountering less matter in the Swiss-cheese model 
than in the EdS cheese model. Moreover, this becomes increasingly true with the 
formation of high-density regions as illustrated in Fig. 14.261 by the evolution of the 
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Figure 4.26: At the top is the evolution of the energy density for the Eds cheese model 
(higher curve) and for the phenomenological model with wq = —1.95 and Wa = 4.28. In 
the middle is the constant-time line averaged density as a function of time for the Swiss- 
cheese model (lower curve) and the cheese-EdS model (higher curve). At the bottom is 
their ratio of the last two quantities as a function of time. 

ratio of the previously calculated average density: it decreases by 17% from the 
starting to the ending time. 

The calculation of Eq. (14.211) is actually, except for some factors like the cross- 
section, the opacity of the Swiss-cheese model. Therefore, a photon propagating 
through the Swiss-cheese model has a different average absorption history; that is, 
the observer looking through the cheese will measure a different flux compared to the 
case with only cheese and no holes. For the moment, in order to explore the physics, 
let us make the approximation that during the entire evolution of the universe, the 
matter is transparent to photons. 

From the plots just shown we can now understand the reason for the best fit 
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values of Wq = —1.95 and Wa = 4.28 found in the case of holes of = 350 Mpc. 
We are using a homogeneous phenomenological model, which has at the present 
time the density of the cheese (see Fig. 14.261) . We want to use it to fit the line- 
averaged density of the Swiss cheese, which is lower than the (volume) averaged 
one. Therefore, going backwards from the present time, the phenomenological model 
must keep its density low, that is, to have a small w. At some point, however, the 
density has to start to increase, otherwise it will not match the line-averaged value 
that keeps increasing: therefore w has to increase toward 0. It is very interesting 
that this simple mechanism mimics the behavior of the concordance- mo del equation 
of state. We stress that this simple mechanism works thanks to the set-up and 
fitting procedure we have chosen; that is, the fact that we matched the cheese-EdS 
solution at the border of the hole, the position of the observer, and the observer 
looking through the holes. Moreover, we did not tune the model to achieve a best 
matching with the concordance model. The results shown are indeed quite natural. 



4.6.6 Beyond spherical symmetry 
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Figure 4.27: Flow chart regarding relationships between the results obtained. 



We have summarized the relationships between the results obtained in this chap- 
ter till now in the flow chart of Fig. 14.271 

Regarding dL{z), we found no important effects from a change in the redshift: 
the effects on diiz) all came from dA driven by the evolution of the inhomogeneities. 



76 



CHAPTER 4. SWISS CHEESE 



Regarding light-cone averages, we found no important effects with respect to the 
expansion: this negative result is due to the compensation in redshift discussed in 
Sec. 14.4. II and it is the same reason why we did not find redshift effects with di{z). 
This is the main limitation of our model and it is due ultimately to the spherical 
symmetry of the model as explained in Sec. 14.4. 1[ 

We found important effects with respect to the density: however this is not due 
to the effects driving the change in dA_- The latter is due to structure evolution while 
the former to the presence of voids, so the two causes are not directly connected. 
Indeed, it is possible to turn off the latter and not the former. 

We can therefore make the point that the expansion is not affected by inhomo- 
geneities because of the compensation due to the spherical symmetry. Density, on 
the other hand, is not affected by spherical symmetry, so there are no compensa- 
tions, and the photon will systematically see more and more voids than structures. 
We can therefore argue that the average of density is more relevant than the average 
of expansion because it is less sensitive to the assumption of spherical symmetry, 
which is one of the limitations of this model. 

The next step is to define a Hubble parameter from this average density: oc 
{p)-y In this way we are moving from a Swiss cheese made of spherically symmetric 
holes to a Swiss cheese without exact spherical symmetry. The correspondence is 
through the light-cone averaged density which, from this point of view, can be seen 
as a tool in performing this step. See Fig. 14.281 for a sketch. 




ip) 




Figure 4.28: Using (p)^ as a tool to go from a Swiss-cheese model with spherically symmet- 
ric holes to a Swiss-cheese model with non-spherically symmetric holes. The insensitivity 
of density to spherical symmetry makes us think that these two models would share the 
same light-cone averaged density. 
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Summarizing again: 

• We started out from a Swiss-cheese model containing spherically symmetric 
holes only. A photon, during its journey through the Swiss cheese, undergoes 
a redshift that is not affected by inhomogeneities. However the photon is 
spending more and more time in the voids than in the structures. The lack of 
an effect is due to spherical symmetry. We focused on this because a photon 
spending most of its time in voids should have a different redshift history than 
a photon propagating in a homogeneous background. 

• Since density is a quantity that is not particularly sensitive to spherical sym- 
metry, we tried to solve the mismatch by focusing on density alone, getting 
from it expansion (and therefore the redshift history). 

• We ended up with a Swiss-cheese model with holes, which are actually not 
spherically symmetric. In this model there is an effect on the redshift history 
of a photon due to the voids. 

• In practice this means that we will use the phenomenological best-fit model 
found, that is, we will use a model that behaves similarly to the concordance 
model. 



4.6.7 Motivations 

Let us go back to the discussion of Sec. 14. 6. 4^ that is, if it is observationally mean- 
ingful to consider light-cone averages of Q as the basis for the correspondence. For 
example, domain averages at constant time are not directly related to observations. 

Here, we are not claiming that light-cone averages are observationally relevantlfj. 
Rather, we are using light-cone averages as tools to understand the model at hand. 
The approach has been explained in the previous section. 



■^However, a density light-cone average is an indicator of the opacity of the universe and, there- 
fore, could be observationally relevant, as explained in the discussion around Fig. 14.261 
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4.7 Renormalization of the matter equation of state 

In this section we will study how the parameters of the phenomenological model 
depend on the size of inhomogeneities, that is, on the size of the hole. We sketched 
in Fig. 14.291 our set-up: we keep the comoving position of the centers of the holes 
fixed. The observer is located in the same piece of cheese. 




Figure 4.29: Sketch of how the size of the inhomogeneity is changed in our model. The 
shading mimics the initial density profile: darker shading implies larger denser. The 
uniform gray is the FRW cheese. The photons pass through the holes as shown by the 
arrows and are revealed by the observer whose comoving position in the cheese does not 
change. The size of the holes correspond to n = 0, 2, 5 of Eq. Ii4.22\) . 

We changed the radius of the hole according to 

= (4.22) 

where Vh is the radius we have been using till now, the one that results in the holes 
touching. The choice of the 1.4 in the scaling is only for convenience. We let n run 
from to 7. 
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In this analysis we will use instead of the energy density in Eq. fl4.17l) . an energy 
density in which only one effective source appears, and the effective source evolves 

as 



^ = (i + ^)3(i+-o«+0 exp (-3w^-±-] with w^'iz) = w^ + w^ 

(4.23) 

We put i? as a superscript on the equation of state in order to differentiate the 
parametrization of Eq. f l4.23p . which we are now using to study renormalization, 
from the parametrization of Eqs. f l4.16M.lTl) . which we used to compare the phe- 
nomenological model to the concordance model. We are not disentangling different 
sources in Eq. f l4.23p because we are interested in the renormalization of the matter 
equation of state of the cheese, that is, on the dependence of upon the size of the 
hole. To this purpose we need only one source in order to keep track of the changes. 

As one can see from Fig. 14.301 we have verified that = for rh — > 0, i.e., we 
recover the EdS model as the best-fit phenomenological model. 

We are interested to see if the equation of state exhibits a power-law behavior 
and, therefore, we use the following functions to fit Wq and w^: 

PO 





= Qo 1 


f Thin) 


<(0) 


VMO) 


Wain) 
<(0) 


= Qa 1 


( Thin) 
U/.(0) 



(4.24) 

We performed a fit with respect to the logarithm of the above quantities, the result 
is shown in Fig. I4.3U[ We found: 

Po = Pa ^ 1-00 

go = g„~ 0.88. (4.25) 
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Figure 4.30: At the top, dependence of Wq (lower points denoted by circles) and 
(upper points denoted by x) with respect of the size of the hole. At the bottom, fit as 
explained in the text. Recall that is today 350 Mpc. 
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Summarizing, we found three important facts. 

• The parameters of the equation of state as a function of the size of the hole 
exhibit a power-law behavior. 

• The power-laws of Wq and have the same scaling exponent. This is actually 
a check: once a physical quantity exhibits a power-law behavior, we expect 
that all its parameters share the same scaling exponent. 

• The scale dependence is linear: the equation of state depends linearly on the 
length of holes the photon propagates through. We stress that the dependence 
we are talking about is not on the scale of the universe, but on the size of the 
holes. 

We can finally ask which size of the holes will give us a phenomenological model 
able to mimic the concordance model. We found that for n = 1, that is for a holes 
of radius = 250 Mpc, we have Wq = —1.4 and = 0.67, which in terms of 
the energy density parametrization of Eq. f l4.17p . corresponds to wq = —1.03 and 
Wa = 2.19. 
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Conclusions 



The aim of this thesis was to understand the role of large-scale non-linear cosmic 
inhomogeneities in the interpretation of observational data. We are stressing the 
word observational because we think that it is the guiding idea to properly set the 
back-reaction problem. 

Back-reaction has indeed many faces and possible interpretations. We split up 
our efforts into two different directions of research. The first is about building a 
realistic (even if still a toy) and exact inhomogeneous universe model in order to 
study directly how inhomogeneities affect observables like the luminosity-distance- 
redshift relation. The other direction is about describing an inhomogeneous universe 
by means of an effective homogeneous phenomenological model, that is, to average 
inhomogeneities. We think that both approaches are worth studying. 

About the latter, we explored two different paths, using two different models re- 
spectively. The first averaging approach is based on Buchert's work. We introduced 
it and discussed it in Chapter [2] and we worked out a model to test it in Chapter [3l 
We think that this framework is very interesting and gives consistent approach to the 
averaging problem. However, the back-reaction source it gives is not clearly related 
to the observed dark energy. It features constant-time averages that do not match 
with an observational approach based on the light cone. Anyway, this scenario gives 
precious hints toward the understanding of the back-reaction problem and allows us 
to check it quantitatively. 

The model we built is an approximate "homogeneous" universe model in which the 
LTB free functions were chosen in order to not single out the center. It was shown 
indeed that it is possible to fit the observed luminosity-distance-redshift relation by 
adjusting the LTB free functions. To achieve this result, however, it is necessary to 
place the observer at the center of a rather big underdensity. Even though we built 
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nothing more than a toy model, we tried to avoid this fine-tuning. 
We found a neghgible back-reaction. We traced its cause back to the spherical sym- 
metry of the model, which is the main limitation of the LTB metric. Indeed spherical 
symmetry introduces self-averaging in the equations, which erases any effect in the 
fiat case where the back-reaction is exactly zero. The curved case then does not de- 
part significantly from the fiat case. We think that the lack of back-reaction effects 
comes from a inefficient use of LTB models and that both the Buchert framework 
and LTB models can give interesting insights. 

In Chapter HI therefore, we built another model, a Swiss-cheese model where the 
cheese consists of the usual Friedmann-Robertson- Walker (FRW) solution and the 
holes are constructed out of a LTB solution. We focused on a Swiss-cheese model 
because, even if it is made of spherical symmetric holes, it is not a spherical sym- 
metric model as a whole. It is a first step to go beyond spherical symmetry which is 
the main limitation of LTB solutions. In order not to fine-tune the position of the 
observer we placed the latter in the cheese and have the observer look through the 
Swiss-cheese holes. 

We first focused on the effects of inhomogeneities on photons. The observables on 
which we focused are the change in redshift Az{X), in angular-diameter distance 
AdA{z), in the luminosity distance-redshift relation dL{z), and in the distance mod- 
ulus Am{z). We found that redshift effects are suppressed when the hole is small 
because of a compensation effect acting on the scale of a hole, due to spherical sym- 
metry, and on the scale of half a hole, due to the matching to the cheese metric. 
The latter is somewhat similar to the screening among positive and negative charges. 
However, we found interesting effects in the calculation of the angular distance: the 
evolution of the inhomogeneities bends the photon path compared to the FRW case. 
Therefore, inhomogeneities will be able (at least partly) to mimic the effects of dark 
energy. We stress that this non-trivial result depends on the presence of a faster- 
than-cheese expanding void which, we think, is the crucial ingredient in studying 
inhomogeneous cosmological models. 

After having analyzed the model from an observational point of view, we set up 
the fitting problem - our second approach to the averaging problem - in order to 
better understand how inhomogeneities renormalize the matter Swiss-cheese model 
allowing us to eschew a primary dark energy. We followed the scheme developed 
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by Ellis and Stoeger, but modified it in order to fit the phenomenological model 
to the Swiss-cheese one. We chose a method that is intermediate between the fit- 
ting approach and the averaging one: we fitted with respect to light-cone averages. 
In particular, we focused on the expansion and the density. While the expansion 
behaved as in the FRW case because of the compensation effect mentioned above, 
we found that the density behaved differently, thanks to its insensitiveness to that 
compensation effect: a photon is spending more and more time in the (big) voids 
than in the (thin) high density structures. This effect is not directly hnked to the 
one giving us an interesting d^. 

The insensitivity to the compensation effect made us think that a Swiss cheese made 
of spherical-symmetric holes and a Swiss cheese without an exact spherical symme- 
try would share the same light-cone averaged density. Knowing the density behavior 
we will, therefore, be able to know the behavior of the Hubble parameter, which, in 
turn, will be the one of the FRW solution with a phenomenological source charac- 
terized by the fit equation of state. In this way we can think to go beyond the main 
limitation of this model, that is, the assumption of spherical symmetry. Prom this 
point of view, light-cone averaged density can be seen as a tool in performing this 
step. 

Then we studied how the equation of state of a phenomenological model with only 
one effective source depends on the size of inhomogeneity. We found that Wq and 
follow a power-law dependence with the same scaling exponent which is equal 
to unity. That is, the equation of state depends linearly on the distance the photon 
travels through voids. 

We finally dealt with the size of the holes in order to mimic the concordance model. 
We found out that the best fit is for holes of radius r/j = 250 Mpc and an equation 
of state with wq — —1.03 and Wa — 2.19. 

We reached the conclusion that it is possible to set the back-reaction problem in 
a physically meaningful way and that it is by no means worth examining the physics 
of the back-reaction in depth. 



CHAPTER 4. CONCLUSIONS 



Appendix A 



Cosmological evolution of Alpha 
driven by a general coupling with 
Quintessence 

We have here presented (see a general model for the cosmological evolution of 
the fine structure constant a driven by a typical Quintessence scenario. We have 
considered a coupling, between the Quintessence scalar and the electromagnetic 
kinetic term F^^^F^^ , given by a general function Bp{(f)). We have studied the depen- 
dence of the cosmological Aa{t) upon the functional form of BF{(f)) and discussed 
the constraints imposed by the data. We have found that different cosmological 
histories for Aa{t) are possible within the available constraints. We have also found 
that Quasar absorption spectra evidence for a time variation of a, if confirmed, is 
not incompatible with Oklo and meteorites limits. 

We have referred to 142] for a general analysis about the fundamental constants 
and their variation induced by a cosmological scalar. 

A.l Introduction 

Over the last few years there has been an increasing interest in the possibility of vary- 
ing the fundamental constants over cosmological time-scales. This has a twofold mo- 
tivation. On one side, several observations point towards the existence of a smooth 
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dark energy component in the universe, which could be modeled via a dynamical 
scalar field called Quintessence (for recent reviews see [SOIESIEII)- In general, we 
expect such a cosmological scalar to couple with some, if not all, the terms in the 
matter-radiation Lagrangian, thus inducing a time variation of physical masses and 
couplings. On the other side, recent improved measurements on possible variations 
of the fundamental constants are opening up the possibility of testing the theoretical 
models to a good degree of precision over a wide range of cosmological epochs. It 
should also be mentioned that, although controversial, some evidence of time vari- 
ation of the fine structure constant a in Quasar absorption spectra was recently 
reported [17] . The cosmological variation of fundamental constants induced by cou- 
plings with the Quintessence scalar is then worth studying in order to see if such a 
field could be responsible for a measurable effect. 

Among all the possibilities, the time-variation of the fine-structure constant is 
the simplest to study both from the theoretical and experimental points of view. 
In this work we have restricted ourselves to this issue. The theoretical study of 
a time- varying fine structure constant dates back to 1982 when Beckenstein [7] 
first considered the possibility of introducing a linear coupling between a scalar 
field and the electromagnetic field. More recently the Beckenstein model has been 
revived, generalized and confronted with updated experimental limits [121 E31 EZ] • 
The concrete case of the Quintessence scalar has been considered too EH [221 
ISU] . However, as we have seen, most authors restrict their studies to the simplest 
case of a linear or quadratic coupling. The possibility of reconstructing the dark 
energy equation of state from a measure of ct-variation has also been proposed in 
the literature [SHISS]- 

In this work we have discussed a general model for the variation of the fine struc- 
ture constant a driven by a typical Quintessence scenario. After briefiy reviewing 
the most recent observational and experimental constraints on the variation of a, 
we have constructed the theoretical framework. In particular we have considered 
the case of a general coupling -Bf(0) (see Eq. OA. 231) below) which includes several 
classes of possible functions. In this way we have been able to study the dependence 
of the cosmological variation of alpha, Aa{t), upon the functional form of Bpiip) 
and discuss the constraints imposed by present data. We have found that different 
cosmological histories for Aa(t) are possible within the available constraints. We 
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have also found that, Webb et al. data ^47^, if confirmed, are not incompatible with 
Oklo and meteorites constraints [2H] • 



A. 2 Overview of the constraints 



Comprehensive reviews about the theoretical and experimental issues connected to 
the time variation of fundamental constants can be found in Refs. [03] , [S2] and [IH] • 
In the following we are summarizing the available constraints on the time variation 
of a, expressed as functions of the redshift z (see also Fig. lA.ip : 

^ a{z) - ctp 
a «o 
where ao = a(0) is the value measured today. 

(1) The most ancient data come from Big Bang Nuclesynthesis (BBN) and give 
[501 E]: 

^ 10^2 z = 10^° - 10*^ . (A.2) 



a 



(2) More recently we have the limit coming from the power spectrum of anisotropics 
in the Cosmic Microwave Background (CMB) [5]: 

Aa 



a 



< 10- 



(A.3) 



(3) From absorption spectra of distant Quasars there are more controversial 
data. Webb and Murphy's groups combined data [17] report a 4 a evidence for a 
variation: Aa/a = (—0.543 ± 0.116) • 10~^ on a cosmological time span between 
z = 0.2 and z = 3.7. This result has not been confirmed by other groups. For 
example, Chand et al. [231 El] give: Aa/a = (-0.06 ± 0.06) ■ 10"^ for z = 2.3 - 0.4 
and Aa/a = (0.15 ± 0.43) ■ 10"^ for z = 2.92 - 1.59. We have chosen to be 
conservative, considering a limit based on the last two results, which is consistent 
with zero variation: 



a 



< 10 



-6 



3-0.4 



(A.4) 



(4) From the analysis of the ratio Re/Os in meteorites dating around 4.56 billion 
years ago it is possible to compute ^^''Re half-life, which gives [52] : 

Aa 



a 



< W 



0.45 



(A.5) 
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(5) From the Oklo natural nuclear reactor that operated 2 billion years ago in 
Gabon, we also have 



Aa 



a 



< 10~' 



z = 0.14 



(A.6) 



(6) We then have limits coming from laboratory measurments which constrain 
the present rate of change of a. Comparing atomic clocks, which use different 
transitions and atoms, what is obtained is [411: 



^ 10"^^ yr"^ z = 
where the dot represents differentiation w.r.t. cosmic time. 



(A.7) 
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Figure A.l: The experimental constraints (l)-(6) discussed above are summarized in the 
picture: log|AQ/a| is plotted as a function of the redshift z. On the right-hand side we 
zoom on z ^ 10. The grey areas are those excluded by present data. 



(7) In addition to the limits discussed above, there is a constraint coming from 
indirect violation of the Weak Equivalence Principle (WEP), parametrized by the 
Eotvos ratio 

rj = 2^1^ (A.8) 
|ai + a2\ 

where ai e 02 are the accelerations of two different test bodies in the Earth gravi- 
tational field. These constraints come from the fact that nucleon masses get elec- 
tromagnetic corrections from quark-quark interactions. As extensively discussed in 
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|33j . the leading term of the electromagnetic contribution comes from the electro- 
static energy of the quark distribution, which is proportional to a. The corrected 
masses, to leading order in a, can then be written as pUl 

rrip = m + a Bp 

rrin = m + a Bn (A. 9) 



where p, n stand for proton and neutron and Bp = 0.63MeV/ cto, Bn = — O.lSMeV/ a^. 
If we suppose that a = a{(p), then we will induce a 0-dependence on the nucleon 
masses: 

6mn = BnSa ; Snip = Bp6a . (A. 10) 



If we define 



J. = ^ = ^B. (A.11) 



we get an indirect violation of the equivalence principle induced by the 'fifth-force' 
mediated by the scalar field 

^ {R^gn + R^gp) (Ai?„^„ + ARpgp) (A.12) 



where: 



Rf = J"' ^ ^ 0.5 Ai?, = '""^'^ ""''^l ^ 0.06 - 0.1 , (A.13) 

+ nn + Up 

and m ~ 931 MeV is the atomic mass unit. The suffix E refers to the Earth, while 
1 and 2 refer to two test bodies having equal mass but different composition. From 
Eqs. flA.lip - flA.i2p we see that any model of a- variation will induce a characteristic 



9p,n 7^ and hence WEP violation: while the first two factors in Eq. flA.12p are 



universal and depend on the Earth composition, the third term is not zero if and 
only if gp^n 7^ and the test bodies have different composition in neutrons and 
protons. The current limits on WEP violations impose [6]: 



7] < 10" 



(A.14) 
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A. 3 The theoretical framework 



Following Olive et al. [53], the most generic action involving a scalar field, the 
Standard Model fields and an hypothetical Dark Matter particle Xy can be written 



as 



S 



1 
4 



\ J d'x^g BF{<j>)F,,F^^ -\j d'x^g 
+ / d^xyf^ ^[%ljjII)4!j + iBj{(t))mj%{)j%ljj\ 

+ j d^x^ [xH-BM^xX^x] (A.15) 

where p= and for the electromagnetic term, for example, = — ^GqA^. 

The index i = 1, 2, 3 refers to the SU (3) gauge group of the Standard Model and j 
runs over the various matter fields. 

The form of the action flA.lSp follows from supplying 0-dependent factors to all 
mass and kinetic terms to the standard Lagrangian (which would have all Bi = 1). 
In general we would expect that all of the i?j(0)'s are switched on, if not forbidden by 
any symmetry principle. However, the theoretical treatment of the full Lagrangian 
is very cumbersome and so the coupling functions Bi{(j)) are usually switched on one 
at a time. In this way one can also disentangle the effects due to each single term. 
Since our focus is on the fine-structure constant a, we will keep only Bp{(f)) ^ 1 and 
set all the other functions equal to 1. 

The relevant part of the action for the effect we have studied is then 



S 



1 



d X\J—g R 



d X\ 
1 

'a 



-9 



d x\ 



-g BFi(j))F^,F^'' 
which allows to define an "effective" fine structure constant 



a{t) 



(A.16) 



(A.17) 



BFim) 

where ao is the value measured today. From ( jA.171) we obtain the relative variation 
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relevant for each cosmological epoch 

Aa_a{t)-ao l-BF{<p{t)) 



(A.18) 



a ao Bp{(f){t)) 

It can immediately be seen that, depending on the cosmological evolution of (f){t) 
and on the functional form of Bp{(f)), the fine structure constant a could in principle 
have had many possible histories during the life-time of the universe. What possi- 
bilities are allowed by a general coupling -Bi?(0) within the available observational 
constraints is then worth studying. 

The relevant equations governing the cosmological evolution in a fiat universe 
are the following 

^(l + 3u^,)Pi (A.19) 

i 

dV 

(f) + 3H4> + — = (A.21) 
a(p 

where i = m, r, (p runs over the matter (including dark matter), radiation and scalar 
components. The relevant equations of state are Wm = for matter, = 1/3 for 
radiation and as defined in Eq. (lA.22p . It is important to note that the evolution 
equation of the Quintessence scalar flA.2ip does not depend on Bp or its derivatives. 



a 


An 


a 


3M2 


H^ = 









This is due to the fact that the statistical average of the term F^^F^^ over a current 
state of the universe is zero. So the only term that drives (j) during the cosmological 
evolution is the potential V^(0). 

Since we are working under the hypothesis that the scalar field in Eq. (1A.21I) 
is the Quintessence scalar, we should also impose the additional constraints coming 
from Quintessence phenomenology. In particular we choose a runaway potential 
which goes to zero as far as the field (f) rolls to infinity, in accordance with the 
observational data. It is also required that the scalar dynamics gives the correct 
value for the equation of state 

^0 = f T l^iti (^-0.7 today) (A.22) 
0^/2 -I- K (0) 

The most general form for the Quintessence potential involves a combination of a 
power-law and exponential terms [IE]. For the purpose of this work, however, we 
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will consider the simplest case of an inverse power-law potential V{(f)) = M^^^^"", 
which gives a late-time attractor equation of state = —2/{n + 2) during matter 
domination pO]. The potential should also be normalized in order to give the correct 
energy density today (p° ~ Vi^cj)) ~ 2/3 p°): this sets the mass scale M. In what 
follows we have chosen n = 1 in the potential in order to have the correct attractor 
equation of state, and so obtain M ~ ■y/2/3 p° Mp. 

We have checked that choosing different Quintessence potentials gives a subdom- 
inant effect on the cosmological variation of a, with respect to changing the coupling 
function In what follows we will then fix 1^(0) = and study the effect 

of different Bf{4')'s. An interesting study, which is complementary to what is done 
here, is that of Ref. [25] where the effect induced by different Quintessence models 
on the cosmological Aa is examined in detail, while keeping the function Bp{(f)) 
fixed. 

In order to be as general as possible we will consider a function BF{(f)) which 
is a combination of different possible behaviors and characterized by a set of four 
parameters that are allowed to vary freely: 

BFi<P) = (1^) ' [1 - C(0 - 00)"] e-(^-*«^ . (A.23) 

This choice is not motivated by a specific theoretical model, but it is rather a working 
tool for obtaining different functional forms of Bp{(j)) and thus cosmological histories 
of a, according to Eq. flA.lSp . We have chosen a combination of functions (power- 
law, polynomial, exponential) that can be switched on and off at will (depending on 
the values of the parameters e, (, t and q), thus giving rise to a variety of possibile 
Bp{(j))^s. In this way we can carry on a unified discussion of a number of different 
models of a variation. 



A. 4 Cosmic evolution of a 



We have numerically solved the cosmological equations ( lA.19l) -( rA.2ip and then plot 



ted the resulting cosmological history of Aa for various classes of functional forms 
of Bpi^cj)), according to Eq. ( lA.lSp . As already mentioned, for illustrative purposes 
we have chosen a scalar potential V = and initial conditions p^^/p^ = 10^° at 
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z = 10^°. Fig. I A. 2 1 shows the corresponding evolution of the energy densities and of 
the scalar equation of state parameter. 




Figure A. 2: Evolution of the energy densities (left) and scalar equation of state (right) 
for a quintessence model with potential V = 1/4) and initial conditions p^/p^ = 10^'^ at 
z = 10^'^. The dot-dashed line represents the energy density of radiation, the dotted line 
the energy density of matter, the green dashed line the energy density of quintessence and 
the red solid line the attractor. All of the energy densities are expressed in units of the 
present critical energy density . 



Linear coupling 

The simplest case is given by the choice e = r = and q = 1 for the parameters in 
Eq. (1X23]) : 

5^(0) = l-C(0-0o) (A.24) 

This case corresponds to the original Beckenstein proposal [7\, which however did 
not supply a potential to the scalar fielcll]. Copeland et al. [25] give a comprehensive 
discussion on various Quintessence models linearly coupled to the electromagnetic 
field, but assuming Webb et al. data [17] to be correct and imposing on Aa{t) to 
agree with that measure. 

^To be precise, Beckenstein actually invoked an exponential coupling, which however is practi- 
cally equivalent to eq. (IA.24p due to the smallness of ^cf). 
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In our case, the resulting Aa, as defined in Eq. (lA.lSp . is plotted in Fig. IA.3[ 
We tried a number of different values for (, in order to verify in which cases all 
the available experimental constraints were simultaneously satisfied. We found that 
they are all respected for ( < 0.6 ■ 10^®. With this choice, the constraints on the 
violation of equivalence principle and the constraints derived from atomic clocks are 
automatically satisfied: 



4 • 10-2^ < 10" 



13 



a 



4 ■ 10^^^ < 10^^^ yr"^ 



(A.25) 
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Figure A. 3: The logarithm of\Aa/a\ is plotted as a function of Log{z + 1) for Bp{(f)) = 
1 - C(0 - -^o) mth C = 0.6 • 10-5 (dotted line), C = 0.6 • 10"^ (solid line) and C = 0.6 • 10"^ 
(dashed line). On the right-hand side we zoom on z ^ 10. Only the curves not overlapping 
the grey areas are phenomenologically viable. 



Polynomial coupling 

A slightly more complicated case is given by the choice e = r = 0, allowing the 
exponent g to be > 1: 

i?^(0) = l-C(0-0o)' (A.26) 

The case of a quadratic couphng (g = 2) was considered in Ref. [3?], but with the 
additional assumption of a proportionality relation between -Bf(0) and V((j)). 

We have found that the data do not impose any upper limit on the exponent q 
and that increasing q makes it possible to reduce the fine-tuning in (. For example. 
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choosing ( = 10~^ the experimental hmits are respected for q = 6 and the constraints 
on the violation of equivalence principle and the constraints derived from atomic 
clocks satisfied by many orders of magnitude. In Fig. IA.4I we plot Log |Aa/Q;| for 
( = 10^^ with g = 3, 6 and 9, as function of red-shift. 



Log of alpha variation 




-Log z+1 -Log z+1 



Figure A. 4: The logarithm of\Aa/a\ is plotted as a function of Log{z + 1) for Bf{(I)) = 
1 ~ C{4' ~ 4'o)'^ with C = 10""^ and q = 3 (dotted line), q = G (solid line) and q = 9 (dashed 
line). On the right-hand side we zoom on z ^ 10. Only the curves not overlapping the 
grey areas are phenomenologically viable. 



As already mentioned, by increasing the exponent q we can do even better. For 
example, with q = 17 the experimental constraints are satisfied even for = 1, 
as illustrated in Fig. IA.5I It should be emphasized that among all the possibilities 
considered in this work, this choice of the parameters appears to be the most natural 
of all. A notable feature is also the fact that the value of Aa is enhanced in the 
past, with respect to the other cases, becoming closer to the observational limits, 
while falling off very steeply in recent times. 

Power— law coupling 

With the choice ( = t = we obtain the following coupling function: 

BF{<j^)=(J-y . (A.27) 

In this case, it is necessary to fine-tune the exponent e in order to satisfy the data, 
due to the smallness of in the early universe. Keeping e of order one would 
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Log of alpha variation 




-Log z+1 -Log z+1 



Figure A. 5: The logarithm of\Aa/a\ is plotted as a function of Log{z + 1) for Bp{(j)) = 
1 ~ C(0 ~ 4'oY with C = 1 and q = 17 (solid line). On the right-hand side we zoom on 
z ^ 10. Note that all the experimental limits are satisfied without any fine-tuning in the 
parameters of the function Bp{(j)). 



violate even the constraints from BBN. We found that the experimental limits are 
respected for |e| < 4 - 10"''. In Fig. IA.6l we plot Log |Aa/a| as a function of red-shift 
for different choices of e. Note that the sign of Aa depends on the sign of e. With 
the choice e = 4 ■ 10"'', the constraints on the violation of equivalence principle and 
the constraints derived from atomic clocks are automatically satisfied: 

rv 

4 ■ 10-" < 10-^^ yr"^ (A.28) 



21 1 n-13 



~ 4 - 10"^' < 10 

Such a small exponent might look quite unnatural, however Eq. ( ]A.27I) for e ^ 1 is 
equivalent to: 



5^(0) = l + eln(^^J . (A.29) 
In this way the fine tuning is moved from the exponent to the coefficient. 



Exponential coupling 

The choice e = C = of the parameters in (1A.23P gives: 



Bf{(I)) = e-"('^-'^°) . (A.30) 

In this case, if r ^ 1 it is not possible to satisfy all the constraints at the same time. 
Depending on the sign of r, the resulting |Aa/a| becomes too large in the early or 
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Log of alpha variation 




-Log z+1 -Log z+1 

Figure A. 6: The logarithm of \Aa/a\ is plotted as a function of Log{z + 1) for Bp{(l)) = 
(^^y with e = 4 • 10-'^ (dotted line), e = 4 • 10"^ (solid line) and e = 4 ■ IQ-^ (dashed 
line). On the right-hand side we zoom on z ^ 10. Only the curves not overlapping the 
grey areas are phenomenologically viable. 



late universe. For r <^ 1, instead, the coupling function becomes equivalent to the 
linear case: 5^7^(0) = e'^^'^-'^") ^ 1 - r(0 - 0o)- 

Linear and power— law coupling combined 

Now let's consider two factors in (1A.23I) with q = I and r = 0. 



5f(</')= (^^J (l-C(</>-</'o)) • (A.31) 

For an arbitrary choice of ( and e, the resulting Aa is similar to the linear coupling 
or power-law coupling case, depending on which factor dominates. It is instead 
interesting to consider the case ( = 'ye in which the two factors are of the same 
order of magnitude. If 7 > 0, the two factors can contribute in an opposite way and 
it is easy to obtain Aa ~ also at some time in the past. For example, with the 
choice e = 2.4 ■ 10~^ and 7 = 2.2 we obtained the behavior plotted in Fig. IA.7 



m 



which we have varied 7 of 10%. For the choice e = 2.4- 10~^, 7 = 2.2 the constraints 
on the violation of equivalence principle and the constraints derived from atomic 
clocks are automatically satisfied: 



20 ^ 1 n-13 



^ ^ 2 ■ 10"^^ < 10 



a 



ao 



9 ■ 10"^^ < 10"^^ yr"^ (A.32) 
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Log of alpha variation 
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Figure A. 7: The logarithm of\Aa/a\ is plotted as a function of Log{z + 1) for Bp{(j)) = 
(i)' (1 - 7 ^ ("^^ - '^o)) witii e = 2.4 • 10-^ 7 = 2.2 (solid line) and 7 = 2.2 ± 10% (dashed 
and dotted respectively). On the right-hand side we zoom on z ^ 10. Only the cmves not 
overlapping the grey areas are phenomenologically viable. 



Power— law and exponential coupling combined 

Since, as already discussed, the exponential coupling function case is equivalent to 
the linear one, this possibility falls within the previous example. 

Polynomial and exponential combined 

Now let's consider e = and q ^ 1. The case g = 1 is not interesting since the 
two factors would be almost equivalent and the behavior corresponding to the linear 
case. Let's choose then, for example, g = 6: 

Bp(<P) = (1 - C (0 - 0o)') e--(*-<^°) (A.33) 

For recent times (z < 1) the exponential coupling e^^^'^^'^°^ ~ 1 — r(0 — 0o) dom- 
inates, while in the past the two terms can be of the same order and, due to q 
being even, cancel at some time. This is shown in Fig. IA.8I The constraints on the 
violation of equivalence principle and the constraints derived from atomic clocks are 
satisfied. 
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Figure A. 8: The logarithm of \Aa/a\ is plotted as a function of Log{z + 1) for Bp{(l)) = 
(1 - C (</> - (t)of) e-^(''^-'^o) with T = 0.6 • lO^^ andC = 2- 10"^ (dotted line), C = 3.2 • 10"^ 
(solid line) and C = 5 • 10^^ (dashed line). On the right-hand side we zoom on z ^ 10. 
Only the curves not overlapping the grey areas are phenomenologically viable. 



A. 5 Summary and Conclusions 

In this work we have carried out a comprehensive study of the cosmological variation 
of the fine-structure constant a induced by the couphng of the electromagnetic field 
with a typical Quintessence scalar. We have considered a variety of functional 
forms for the coupling function Bp{(f)), obtainable from a general expression (see 
Eq. (IA.23P ) depending on four parameters. 

We have found that very different cosmological histories for Aa are possible, 
depending on which parameters are switched on. For example, we can produce a 
Aa which is well below the observational constraints in the early universe and just 
within the experimental limits in recent times (linear coupling case). But also the 
converse is possible, if we choose a polynomial coupling. In particular, the behavior 
at small redshift can be qualitatively very different depending on the model we 
choose: sharply decreasing in the polynomial coupling case or mildly decreasing 
with the power-law coupling. 

By combining different functional forms, a notable feature emerged. In some 
cases it is possible that the scalar dynamics drives Aa to a zero at some time in the 
past, thus inverting the slope of its cosmological evolution. This happens for all the 
combined cases discussed here. 
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It is also worth remarking that in our parameter space span we have found 
solutions with extremely reduced fine-tuning, which are still compatible with the 
available constraints. This is the case of the polynomial coupling with exponent 
g > 15 lifting the fine-tuning of the coefficient ( (usually constrained to be < 10~®) 
to order 1. 

It should also be emphasized that, while in the literature the result by Webb et 
al. [17] is usually said to be incompatible with the Oklo limit [28], we have found 
that this is not always the case. For example, in the polynomial coupling case it is 
possible to obtain several examples with a Aa which matches the Quasars data and 
at the same time respects the Oklo bound. 



Appendix B 
LTB models 



We will introduce here the Lemaitre-Tolman-Bondi (LTB) solution, which will be 
used in Chapters [3] and HI 

First, in Table IB.ll we list the units we will use for mass density, time, the 
radial coordinate, the expansion rate, and two quantities, Y{r,t) and W{r), that 
will appear in the metric. 



Table B.l: Units for various quantities. We use geometrical units, c = G = 1. Here, 
the present critical density is pco = S-ffg obs/'^'^' with H^^obs = 70 km Mpc~^ . In 
order to have the proper distance today we have to multiply the comoving distance by 
a{to) ~ 2.92. 

Quantity Notation Unit Value 

mass density p{r,t), p{r,t) pco 9.2 x 10~^° g cm~^ 

time t, T, t, tsB, Tq (evrpco)""^^^ 9.3 Gyr 

comoving radial coordinate r (STrpco)"^^^ 2857 Mpc 

metric quantity (Svrpc'o)"^''^ 2857 Mpc 

expansion rate H{r,t) {GTrpcoY^^ |-f^o,06s 

spatial curvature term W{r) 1 — 



The time t appearing in Table ET] is not the usual time in FRW models. Rather, 
t = T — To, where T is the usual cosmological time and To = 2Hq^ /3 is, for an 
EdS model, the present age of the universe. Thus, t = is the present time and 
t = tBB = —To is the time of the big bang. Finally, the initial time of the LTB 
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evolution is defined as t. 

Both the FRW and the LTB metrics, after solving Einstein's equations, can be 
written in the form: 

rfs2 ^ _^^2 ^ T^^dr'' + Y^r, t) dn^ (B.l) 

where the "prime" superscript will denote d/dr and the "dot" superscript d/dt. 
Y{r, t) is the position of the shell r at the time t. In order not to have shell crossing 
we will demand F' > 0. 

It is clear that the Robertson-Walker metric is recovered with the substitution 
y(r, t) = a{t)r and W'^{r) = 1 — kr"^. The above metric is expressed in the syn- 
chronous and comoving gauge. 



B.l The EdS solution 

To connect the LTB solution to the FRW one, let's see how a flat matter-dominated 
universe (the EdS model) is described within the LTB notation. For the EdS solution 
there is no r dependence to p or H. Furthermore, Y{r,t) factors into a function of 
t multiplying r {Y{r,t) = a{t)r), and in the EdS model W{r) = 1. In this model 
Qm = 1, so the value of p today, denoted as po, is unity in the units of Table [BTTI 
In order to connect with the LTB solution, we can express the line element in the 
form: 

els'" = -dt" + Y'\r, t)dr^ + Y''{r, t) d^f (B.2) 

The Friedman equation and its solution are (recall t = corresponds to the present 
time) : 

H\t) = ^p(t) = ^(t + l)-^ (B.3) 

it + 1)2/3 

F(r,t) = ra{t)=r^j^^^ (B.4) 

where the scale factor is normalized so that at the beginning of the LTB evolution 
it is a{t) = 1. 

For the EdS model, Tq = 1. We also note that the comoving distance traveled 
by a photon since the big bang is tbb = ^/do- 
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B.2 The LTB solution 

The LTB model [381 [611 [TO] is based on the assumptions that the system is spheri- 
cally symmetric with purely radial motion and the motion is geodesic without shell 
crossing (otherwise we could not neglect the pressure). 

It is useful to define a "Euclidean" mass M(r) and an "average" mass density 
p(r, t), defined as: 

r An 
M{r) =An p{r, t) Y^Y' dr = —Y^{r, t) p{r, t) (B.5) 
Jo 3 

In spherically symmetric models, in general there are two expansion rates: an an- 
gular expansion rate, = y(r, t)/y(r, t), and a radial expansion rate, Hr = 
Y'{r,t)/Y'{r,t). (Of course in the FRW model Hr = H±.) The angular expansion 
rate is given by: 

Hlir,t) = -p{r,t) + ^;^ (B.6) 
Unless specified otherwise, we will identify H± = H. The integral of flB.6p is: 



"Y / '?r'\ ^^^^ 

t = t{r)+ I (w^-l + — ) dy (B.7) 



'y \ y J 

From the equation of motion flB.6p we see that the spherical shells of matter further 



away from the origin in relation to a particle P, do not affect the motion of P at all. 
This can be traced back to the Birkhoff's theorem. 
Differentiating flB.6p gives the acceleration equation: 

^ = -^p('^,^) (B.8) 



Equations (IB. 61) and (IB.SP are similar to the homogeneous ones with the crucial 



difference that there is an extra dependence on r. 

Here we find the first potential problem of this metric as far as averages are 
concerned: in the equations ( lB.6IIB.8j) the density is already averaged. This is the 
effect of the spherical symmetry. Moreover this "automatic" mean is taken from 
the center of symmetry which remains singled out in spite of our choices of the 
arbitrary functions. This mean that it could be inappropriate to average again the 
LTB model. 
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To specify the model we have to specify initial conditions, i.e., the position 
Y{r,t), the velocity Y{r,t) and the density p(t) of each shell r at time t. In the 
absence of shell crossing it is possible to give the initial conditions at different times 
for different shells r: let us call this time t{r). The initial conditions fix the arbitrary 
curvature function W{r): 



W\r) - 1 ^ 2E{r) = ( _ 1 ^ 



:b.9) 



where we can choose Y{r^t) = r so that M(r) = A-n /J' p(f, t) dr. 

In a general LTB model there are therefore three arbitrary functions: p{r,t), 
W{r) and t(r). In Sect. lB.3l we provide a discussion about the number of independent 
arbitrary functions in a LTB model. 

Finally here is the expression of the full Ricci scalar R: 



VK^const 

F2 Y Y' YY' - 1 ^~WW 



FLRW LTB VK^l 

and of the spatial Ricci scalar TZ: 

W'^ - I WW E E' 

^ = -2^7^ ^YY^ = -^Y-^-^YY' 

and of the expansion rate: 

^ = = 2^ + ^ (B.12) 

and of the square shear: 



k 
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B.3 About the arbitrary functions in a LTB model 

Here we illustrate, by means of an example, the choice of the arbitrary functions in 
LTB models. We are going to analyze the flat case. Indeed we have an analytical 
solution for it and this will help in understanding the issues. 

We said previously that there are three arbitrary functions in the LTB model: 
p(r), W{r) and t(r). They specify the position and velocities of the shells at a 
chosen time. In general, t depends on r; because of the absence of shell crossing it 
is possible to give the initial conditions at different times for different shells labeled 
by r. 

We start, therefore, by choosing the curvature E{r) = {W'^{r) — l)/2 to vanish, 
which can be thought as a choice of initial velocities Y at the time t(r): 



(B.14) 



r, t(r) 



For E{r) = 0, the model becomes 

ds'^ = -dt^ + dY^ + Y^dn\ (B.15) 



with solution 



p{r,t) = [^-^(rr^ (B.16) 



where 



p(r,t(r)) 



t(r) = t{r) - p'^^^{r,t{r)) 
3M(r) 1 



An F3 



(B.17) 



r, t(r) 

The next step is to choose the position of the shells, that is, to choose the density 
profile. As far as t(r) is concerned, only the combination i{r) matters. This, however, 
is not true for M(r), which appears also by itself in Eq. ( ]B.16|) . 



Looking at Eq. ( ]B.17p we see that to achieve an inhomogeneous profile we can 
either assign a homogeneous profile at an inhomogeneous initial time, or an inho- 
mogeneous density profile at a homogeneous initial time, or both. Moreover, if we 
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assign the function M(r), then we can use our freedom to relabel r in order to ob- 
tain all the possible i{r). So we can see that one of the three arbitrary functions 
expresses the gauge freedom. 

In this thesis we fixed this freedom by choosing i{r) = t and Y{r,t) = r in order 
to have a better intuitive understanding of the initial conditions. 



Appendix C 

Sewing the photon path 



In this Appendix we will show how to sew together the photon path between two 
holes. We will always use center-of-symmetry coordinates, and therefore we will 
move from the coordinates of Oi to the ones of O2 illustrated in Fig. IC.ll The 
geodesic near the contact point G is represented by the dashed line segment in Fig. 





Y 













Figure C.l: Illustration of the procedure to calculate the transition between two 
holes. The dashed line is a segment of the geodesic. Oi and O2 represent the two 
coordinate systems. 

First, we want to find out the value A of the affine parameter for which the 
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photon is at G. This is found by solving 

G2 = (ri(A) cos0i(A) - 2rft, ri(A) sin0i(A)) 
rl = xIg + vIg- (C.l) 

These equations imply 

rliX) + 3rl - 4r,{X)rh cos0i(A) = 0. (C.2) 

Then we can give the initial conditions for the second hole: 

92(A) = gi(A) 

t2(A) = ti(A) 

'^2(A) = rh 

02(A) = arccos(x2G/^/i)- (C.3) 

Finally, we will need the constant c^, a sort of constant angular momentum density. 
Repeating the procedure of Sect. 14.3.11 for the first hole we will find 

C20 = sina2 gi(A) FsIa • (C.4) 

Only 0^2 is missing. One way to find it is to calculate the inner product in Oi 
coordinates of the geodesic with the normalized spatial vector parallel to O2G (see 

Fig.m]). 
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